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ABSTRACT 


This work is a logical continuation of an earlier study (NAS 8-32614) 
on electrophoresis technology and forms part of a two-pronged effort 
involving research at Princeton University and the Marshall Space Flight 
Center (MSFC) . The overall objective is to establish theoretically and 
confirm experimentally the ultimate capabilities of continuous flow 
electrophoresis chambers operating in an environment essentially free of 
particle sedimentation and buoyancy. The efforts at Princeton are devoted 
to : 


(i) Studying the effects of particle concentration on sample 
conductivity and dielectric constant. The dielectric constant 
and conductivity have been identified as playing crucial 
roles in the behavior of the sample and, thus, on the 
resolving power and throughput of continuous flow devices. 

(ii) Improving the extant mathematical models to 
predict flow fields and particle trajectories in continuous 
flow electrophoresis. 

A dielectric spectrometer was designed and built to measure the complex 
dielectric constant of a colloidal dispersion as a function of frequency 
between 500 Hz and 200 kHz. The real part of the signal can be related to 
the sample's conductivity and the imaginary part to its dielectric constant. 
Measurements of the dielectric constants of several different dispersions 
disclosed that the dielectric constants of dilute systems of the sort 
encountered in particle electrophoresis are much larger (e n 1000 ) than 
would be expected based on the extant theory. Experiments were carried out 
to show that, in many cases, this behavior is due to the presence of a 
filamentary structure of small "hairs" on the particle surface. 

A technique for producing "electrokinetically ideal" synthetic latex 
particles by heat-treating was developed. Here the term "electrokinetically 
ideal" means that the particles conform to the tenets of the classical 
theory. 

Given the ubiquitous nature of "hairy" surfaces with both cells and 
synthetic particles, it was deemed necessary to develop a theory to explain 
their behavior. A theory for electrophoretic mobility of hairy particles 
was developed; an extension to include conductivity and dielectric behavior 
is underway. 

Finally, the extant computer programs for predicting the structure of 
electro -osmotically driven flows were extended to encompass flow channels 
with variable wall mobilities. 




SUMMARY AND CONCLUSIONS 


Earlier work at Princeton (under contract NAS 8-32614) showed that 
particle concentration has a strong influence on sample conductivity (more 
specifically, on the "real" or "DC" conductivity) . Concurrently, 
experimental studies by P. H. Rhodes at MSFC showed that sample and buffer 
conductivities profoundly influence spreading of the sample. A large 
conductivity mis -match between sample and buffer causes sample to spread 
rapidly from the front to rear walls of the channel. Rhodes developed an 
electrohydrodynamic theory of this spreading which shows that in addition to 
the conductivity, the dielectric constant should also affect spreading 
behavior. Therefore, to optimize performance of a continuous flow device it 
will be necessary to understand the spreading process and how it is 
influenced by the conductivity and dielectric constant of the sample. 
Accordingly, one needs to be able to measure the dielectric behavior of 
particulate samples as well as their DC conductivity. 

Much of the work accomplished under the present contract has been 
published in the refereed literature; copies of the papers that have 
appeared or are "in the press" are attached and copies of the two PhD 
theses are available from the principal investigator. Accordingly, only a 
brief synopsis of the results is given here. 

The dielectric spectrometer designed and built at Princeton measures 
the properties of a given suspension at frequencies between 500 Hz and 200 
KHz^-»^. A complete frequency scan takes about 25 minutes. The instrument 
has been used to measure the dielectric constant and conductivity of latex 
suspensions' 1 » J , similar to those used in Rhodes' experiments, and colloidal 
silica particles^. The results are linear in the volume fraction of 
particles for dilute systems (< 8 vol%) as they should be. One surprising 
result is that the low frequency limit of the dielectric constant is much 
larger (over 1000 for some suspensions) than the extant theory predicts. In 
addition, the relaxation frequency is lower than predicted by the theory. 

It follows from these results that the effect of the dielectric 
constant on the sample's configuration in an electrophoresis machine may be 
even larger than anticipated. Clearly, it is of crucial importance to 
understand these phenomena if electrophoresis is to be used to fractionate 
cells or cell fragment populations. 

Inasmuch as the dielectric constants of a number of suspensions were 
found to exceed those predicted by the classical theory, it is natural to 
inquire as to why. According to the classical theory the particle surfaces 
are smooth and the charge is immobile. Thus, one ought to be able to 
measure the particle mobility, deduce the charge from the classical theory, 
and use it to predict the dielectric response. This turns out not to be the 
case is most situations. The weak point in the theory appears to be in the 
notion of an immobile charge. For example, if the interfacial region is 
somewhat thick, as it would be if the shear surface is displaced a small 
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Numerals refer to the papers in the bibliography. 
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distance from the particle surface, then ions might migrate laterally within 
this region. Moreover, it is believed that many synthetic latex particles 
are hairy due to the presence of (i) oligomeric chains adsorbed or otherwise 
attached to the surface or (ii) polymer filaments protruding from the 
surface. Hairiness provides a means for moving the shear plane out away 
from the particle and furnishes a relatively thick region in which ions 
might move across the surface in response to lateral forces. 

We sought to test the hairy layer hypothesis by suppressing or removing 
the layer. This was done with two sorts of particles: native surfaces with 
naturally occurring hairiness and with particles to which covalently bound 
polymer filaments are attached^ ’ . A mild heat treatment was used wherein 
particles were heated to 120°C for 1 to 12 hours. We were able to follow 
the collapse of the hairy layer using measurements of electrophoretic 
mobility and dielectric behavior. It should be noted that the dielectric 
measurements were much more sensitive to changes in the surface morphology 
than measurements of particle mobility. With some particles, our annealing 
procedure produced "electrokinetically ideal” particles where the dielectric 
response could be predicted accurately from the surface charge inferred 
from the mobility measurements. Moreover, the dielectric constants of such 
suspensions are an order of magnitude smaller than those comprised of native 
particles. It would be interesting to carry out some of Rhodes' 
electrohydrodynamic experiments with such particles . 

As matters now stand, it appears that the dielectric constants and 
conductivities of most synthetic particle dispersions are heavily 
influenced by the presence of a hairy surface texture. This may also be the 
case with many biological cells and their dielectric properties should be 
measured with the new spectrometer. Clearly, surface texture affects the 
electrokinetic response of any dispersion. 

There is, at present, no theory which encompasses the behavior noted 
above. As part of the current effort, a theory for the electrophoretic 
mobility was developed^ > ® . Work on the corresponding theory for dielectric 
behavior is underway^ . 

Under the auspices of the previous contract some work was done on the 
mobilities of particles in concentrated suspensions. Two papers were 
completed during the present contract dealing with that subject^ 1 

Finally, we also updated our program for predicting the flow structure 
in continuous flow devices to encompass non-uniform electro -osmosis along 
the boundary of the flow channel . 
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Dielectric Spectroscopy of Colloidal Suspensions 

I. The Dielectric Spectrometer 

D. F. MYERS AND D. A. SAVILLE 1 
Department of Chemical Engineering. Princeton University, Princeton. New Jersey 08544 
Received July 11, 1988; accepted November 7, 1988 

Four-electrode signal-processing techniques have a number of advantages over conventional bridge 
techniques in dielectric spectroscopy, but are not widely used because complex electronic circuits are 
required. A new design for a dielectric spectrometer of this type which circumvents these problems is 
described here. The device has been tested using electrical networks, pure electrolytes, and dilute colloidal 
suspensions. In experiments with a dilute, amphoteric latex, large suspension dielectric constants are 

found. C 1989 Aademic Press, Inc. 


I. INTRODUCTION 

Dielectric spectroscopy offers certain ad- 
vantages over the static techniques customarily 
used to study the electrokinetic properties of 
colloidal particles. Measuring the dielectric 
constant and conductivity over a range of fre- 
quencies provides information on temporal 
processes in the Stem and diffuse layers as- 
sociated with the particles. In addition, the low 
frequency limits of the measurement can be 
related to static measurements. Unfortunately, 
it is difficult to measure the dielectric behavior 
of dilute, conducting suspensions, and no “off- 
the-shelf” instrument for this purpose exists. 
The most serious experimental problem is 
electrode polarization, which arises whenever 
current passes through an electrode-electrolyte 
interface. In the simplest measurement con- 
figuration, current passes through a sample 
placed between two electrodes and the voltage 
drop divided by the applied current yields the 
complex impedance. However, electrode po- 
larization creates an additional voltage drop 
at the electrode-electrolyte interface which is 
difficult to measure or predict. This clouds in- 
terpretation of the results. 

1 To whom correspondence should be addressed. 


Another problem stems from the need to 
study dilute suspensions so that the experi- 
mental data can be interpreted using rigorous 
theories. Since the contribution of the particles 
to the suspension ’s electrical properties is pro- 
portional to the volume fraction of particles, 
the signal-to-noise ratio is low in dilute sus- 
pensions. Also, small temperature fluctuations 
in the electrolyte produce changes in the elec- 
trolyte conductivity larger than the change due 
to the particles, creating a stringent tempera- 
ture control requirement. 

Nevertheless, dielectric spectroscopy has 
been used successfully with colloidal disper- 
sions. Schwan el al. ( 1 ) used apparatus con- 
sisting of an AC Wheatstone bridge and a 
variable electrode-gap cell, the latter to mini- 
mize errors due to electrode polarization. Di- 
electric constants as large as 2500 were ob- 
served with 30 vol% suspensions of 0.188-/xm 
polystyrene particles. Springer et al. (2) and 
Lim and Franses (3) employed cells similar 
to those of Schwan et al., with a transformer 
ratio-arm bridge in place of the Wheatstone 
bridge. Large dielectric constants were mea- 
sured although the suspensions were much 
more dilute than those used by Schwan et al. 
Recently, Midmore et al. (4) used a system 
much like Schwan ’s to measure electrical 
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conductivities at low frequency. Midmore and 
O'Brien (5) measured high-frequency dielec- 
tric spectra of concentrated suspensions with 
a fixed-electrode cell and a bridge circuit. 

Despite its record of success. Schwan’s two- 
electrode technique (6) has a number of 
drawbacks. First, the method requires several 
measurements at every frequency, each re- 
quiring manual adjustment of the electrode 
spacing. This limits the degree to which the 
measurement can be automated and increases 
the measurement time substantially. Second, 
the signal-to-noise ratio is poor when dilute 
suspensions are involved. To eliminate errors 
due to electrode polarization and avoid man- 
ual adjustments or multiple measurements at 
each frequency. Hayakawa et al. (7) used a 
four-electrode cell design originally proposed 
by Schwan (6). An analog and digital signal- 
processing scheme improved the signal-to- 
noise ratio of the measurement and an elec- 
tronic compensation circuit removed the need 
for precise temperature control. How-ever. the 
technique is difficult to implement due to the 


complexity of the circuitry and the need for 
very precise electronic components. 

II. THE NEW INSTRUMENT 

We set out to develop a device which would 
incorporate the advantages of the Hayakawa 
instrument (7) in a simpler design. A four- 
electrode cell geometry was chosen, w'ith a de- 
tector patterned after those used by Hayakawa 
el al. (1) and Nakamura el al. ( 8 ) . The signal 
generator and signal processing circuits were 
based on a new design which allows the use 
of off-the-shelf equipment and facilitates 
computerized control and data acquisition. 

Our instrument employs a four-electrode 
sample cell (Fig. 1 ). The total volume of the 
sample contained in the cell is approximately 
2 cm 3 . In Fig. 2, current (7) enters the cylin- 
drical cell through the source electrode and 
exits at the sink electrode, each having large, 
disk-shaped surfaces. Two needle-shaped 
sensing electrodes measure the voltage drop 
across the middle of the cell. Platinum foil was 



Fig. 1. A photograph of the four-electrode sample cell. 

Journal of Colloid and Interface Science. Vol. 1 3 1 . No. 2. September 1 989 
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Fig . 2. A diagram of the four-electrode cell configura- 
tion. 


used for the source-sink electrodes and plat- 
inum wire for the sense electrodes. All elec- 
trodes were platinum blacked using Schwan’s 
procedure (6). Electrode polarization influ- 
ences the sense electrode voltage only to the 
extent that the voltage measurement causes a 
current flow. Current flow is minimized by 
connecting each sense electrode to a high- 
impedance, unity-gain buffer amplifier. A dif- 
ferential amplifier subtracts the “negative” 
sense voltage Vi from the “positive” sense 
voltage V\ to yield the voltage drop A V. Thus, 
if the magnitude of the current (/) passing 
through the cell via the source electrode is 
known, the impedance in the region between 
the two sense electrodes can be computed by 
simply dividing the voltage drop by the cur- 
rent, Z = AF/7. Alternatively, the admittance 
is given by Y - 7/ AF. 

This technique is insensitive to the effects 
of electrode polarization because the four- 
electrode configuration separates electrode 
functions. No single electrode is required to 
pass current and measure voltage simulta- 
neously. Electrode polarization at the source 
and sink electrodes affects the voltage drop 
across the entire cell, but not the current, 
keeping 7 unchanged. Although electrode po- 
larization at the sense electrodes is large, the 
high impedance of each buffer amplifier min- 
imizes current flow through the electrodes so 
the voltage drop is negligible compared to AF. 

Journal of Colloid and Interfact Science, Vol. 131, No. 2, September 1989 


The complete instrument consists of the 
four-electrode cell, signal generator, detector, 
and automatic balancing circuit sections, along 
with two lock-in amplifiers (Fig. 3 ). The signal 
generator section provides an alternating cur- 
rent to the source electrode using sine-wave 
generators, a summing amplifier, and a volt- 
age-to-current (F/7) converter. Two sine 
waves are produced by the generators, one at 
a measuring frequency a> and the second at 
the reference frequency u> r . The summing am- 
plifier combines the two signals to produce a 
composite waveform, and the V/l converter 
circuit changes the constant-amplitude voltage 
signal to a constant current 

/ = 7(o>) + /(<*>,) [1] 

independent of the impedance presented by 
the cell. This ensures that the current in the 
cell and A V will be independent of fluctuations 
in the source electrode impedances. 

Buffer amplifiers, a differential amplifier, a 
summing amplifier, and a current-to-voltage 
converter form the detector. It was noted ear- 
lier that simply measuring the differential 
voltage and the cell current is sufficient for a 
determination of the complex admittance. 
However, instead of measuring the current 
through the cell directly, it is useful to produce 
a voltage proportional to the current. This is 
the function of the current-to-voltage (7/F) 
converter. The output of this circuit is 

F f = -7Z f , [2] 

where Z F is the impedance of the feedback 
network consisting of a resistor R ? in parallel 
with a capacitor C F . Note that the “+” and 
terminals of the op-amp draw no current, 
so the current passing through the feedback 
network (7? F , C F ) is 7, the current which passes 
through the sample. 

The summing amplifier combines the out- 
puts of the 7/F converter and the differential 
amplifier to form the inverted sum 

V 0 - — (AK + Kf), [3) 

representing the balance voltage. If Zq rep- 
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signal generator section 



Fig. 3. Block diagram of the four-electrode frequency-difference dielectric spectrometer. 


resents the impedance of the region of the cell 
between the sense electrodes, then 

AV = IZc [4] 

and 

V 0 = -nZc-Zy). [5] 

This is the final output from the device: a volt- 
age proportional to the difference between the 
sample impedance and the impedance of the 
reference network. A lock-in amplifier mea- 
sures the in-phase and out-of-phase compo- 
nents of this signal. 

The reason this approach increases the sig- 
nal-to-noise ratio of the experiment can be 
seen from the following argument. Recall, first, 
that the signal applied to the celi is composed 
of sine waves at the reference frequency o> r and 
the measuring frequency u>. Suppose, now, that 
the measuring lock-in amplifier is tuned to a? r . 
Then, R F and C F can be adjusted until both 
the in-phase and the out-of-phase components 
of the output signal at u> T are nulled. From Eq. 


[ 5 ) we see that a null output at o> r means that 
Z c (<D r ) is equal to Z F (w r ). 

Next note that the admittance of a material 
in a cell with cell constant G is 

Y = G(a + /wee o)- [6] 

In the case of a pure electrolyte, a and c can 
be considered frequency-independent con- 
stants, a* and c* at the frequencies employed 
here; to is the permittivity of the vacuum. 
Thus, the admittance 

Y = G(ffec + iWEooto) [7] 

has a real part which is constant, and an imag- 
inary part linear in the frequency o>. In a sim- 
ilar manner, the admittance of an electrical 
network consisting of a resistor and a capacitor 
in parallel can be written 

y = i + <«c, [81 

where R is the resistance and C is the capac- 
itance. Note the similarity between the RC 
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parallel circuit admittance and the electrolyte 
admittance. The RC parallel circuit models 
the behavior of the pure electrolyte if com- 
ponent values are chosen such that 



[9] 

C C (. yr (Q , 

[10] 


At null, Eqs. [9] and [10] yield C = C F and 
R = R F . 

Now consider the output when the lock-in 
amplifier is tuned to the measuring frequency 
o) under two experimental conditions: first, 
when the sample consists of a pure electrolyte, 
and, second, when the sample is a suspension. 
With pure electrolyte the sample impedance 
can be accurately modeled as a resistor and 
capacitor in parallel, as shown above. Since 
the reference network is also an RC parallel 
network, the null condition Z c (w r ) = Z F (u r ) 
implies that R F equals the equivalent resistance 
of the sample, and C F equals the equivalent 
capacitance of the sample. Because the fre- 
quency dependences of the electrolyte and the 
RC circuit are the same, Z F (w) = Z c (o>), and 
the output will remain zero at all frequencies. 
The reference network models the frequency 
response of the sample. When the model is 
accurate, as it is when the sample is a homo- 
geneous electrolyte, the output remains null 
at frequencies below 1 GHz. 

With suspension in the cell Z F (cu) will differ 
from Z C M and the output signal will be 
nonzero. The amplitude and phase shift of the 
output reflect the presence of the particles in 
the suspension and will be proportional to 
their volume fraction. Accordingly, the influ- 
ence of the impedance of the background 
electrolyte has been eliminated and a signal 
sensitive to the properties of the particles ob- 
tained. 

The four-electrode cell supresses electrode 
polarization effects and the detector increases 
the signal-to-noise ratio. A third problem 
which must be addressed is the error caused 
by small temperature changes during a fre- 


quency scan. Assume, for the moment, that 
the cell is filled with a pure electrolyte. To 
maintain zero output voltage, the impedance 
Z F must mimic the electrical behavior of the 
pure electrolyte. Temperature changes in the 
electrolyte cause the resistance of the sample. 

R c , to change so R F will not match the sample 
resistance. A mismatch between R F and Rc 
produces an output proportional to the differ- 
ence and reduces the resolution of the instru- f 
ment. If the resistance R F could track changes 
in the electrolyte conductivity, the system 
output would remain null despite thermal drift 
and the error would be eliminated. Tracking 
is performed by the automatic balancing cir- 
cuit (Fig. 4). A second lock-in amplifier, the 
“reference” lock-in, is tuned to the reference 
frequency <u r and monitors system output. The 
lock-in output is a dc signal, proportional to 
the ac amplitude of the in-phase output voltage 
at a> r . This amplitude reflects the amount of 
mismatch between R F and the resistance of 
the sample at the reference frequency and can 
be used to decrease the mismatch by compen- 
sating R f for changes in Rc. The dc output 
from the lock-in is applied to an interface cir- 
cuit which controls the current through a light- 
emitting diode (LED). The LED is placed in 
close proximity to a cadmium sulfide (CdS) 
photoconductive cell, a device which has a re- 
sistance dependent upon the light intensity at 
its surface. Together, the LED and the CdS 
cell form a voltage controlled resistor. The 
LED/CdS combination is placed in series with 
a conventional resistor to form the feedback 



Fig. 4. The automatic balancing circuit. 
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resistance R F . R F can then be varied auto- 
matically by controlling the current flow to 
the LED. 

Suppose Z c (o> r ) = Z F (wy) at the outset of 
an experiment, and that an increase in the 
sample temperature has caused Z c (w r ) to de- 
crease. Whenever the in-phase output becomes 
nonzero, the current through the LED is 
changed so that the difference between R F and 
Rc decreases in magnitude to renull the in- 
phase output. Here, since F 0 (w r ) is positive, 
current flow through the LED increases. This 
reduces the resistance of the CdS cell, thereby 
reducing R F until it is again equal to 7? c (w r ), 
and F 0 (u> r ) vanishes. In this way R F is com- 
pensated for the change in the sample resis- 
tance. Typically, the “effective” level of tem- 
perature variation is less than 0.000 1 °C during 
the 20- to 30-min period required to obtain a 
spectrum when this circuit is used. 

The instrument described here is different 
from that designed by Hayakawa et al. ( 7 ) in 
some important ways. First, custom circuits 
in the signal generator and demodulator sec- 
tions of Hayakawa’s device were replaced by 
commercially available equipment in order to 
facilitate implementation of the technique. 
However, a simple substitution of standard 
equipment for custom made components was 
not possible, so significant design changes were 
required. 

In Hayakawa’s system, the digital signals 
used to create sine waves in the signal gener- 
ator section function as syncronization signals 
in the demodulators, thereby defining the 
point of zero phase shift in an unambiguous 
way. However, phase shifts in the bandpass 
filters of the signal generator create very large 
phase errors which are eliminated by manually 
adjusting the phase characteristics of the filters. 
In our system the function of the demodula- 
tors is performed by the Princeton Applied 
Research (PAR) 124 lock-in amplifier. Instead 
of relying on a fixed phase reference, the phase 
of the system output F 0 is measured relative 
to the phase of AF, the differential amplifier 
output. This is a convenient choice since the 


ratio 1 ' 0 /AFisthe quantity used in data anal- 
ysis, and no manual nulling of phase errors is 
required. 

A second difference between our approach 
and Hayakawa’s is the addition of the balance 
capacitance C F in the detector circuit. The ad- 
dition of this capacitance stabilizes the J/V 
converter circuit without changing the fre- 
quency response of the operational amplifier 
used in the circuit. Also, C F can be adjusted 
so that the in-phase and out-of-phase com- 
ponents of V 0 are similar in magnitude, thus 
minimizing errors which would occur if the 
phase shift of the signal was near 0° or 90°. 


III. DATA ANALYSIS 


The output voltage of the experiment is re- 
lated to the electrical properties of the sample 
through Eq. [5] but the electrical properties 
of the sample are buried in the complex quan- 
tity Zc so V 0 , /, and Z F must be measured. 
Since the sample current (7) is difficult to 
measure directly, it is useful to measure a volt- 
age proportional to it. A particularly conve- 
nient voltage is AV, which was shown to be 
related to 1 through Eq. [4]. Combining this 
with Eq. [5] we obtain 



V 0 and A Fare both measurable quantities. 

The complex impedances Z F and Z c can 
be written in terms of real quantities, 


1 1 + ■ r 

— — + IQ)C F 

z F r f 


1 


= G(a 4- 


[ 12 ] 


[13] 


where G is the cell constant. Substituting these 
expressions into the equations above yields 


a + /ujtfo = 


1 + 


AF 

I 

X 


l*F 


+ /o>cy 


u 


[14] 
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Equating the real and imaginary parts of the 
equation shows that 


ff = 5 Re 


1 + 


AV 


1 

L*f 


+ iojCf 


115] 


and 


€ = — — lm 
Owco 


1 + 


±L 

AV 


L*f 


+ I^Cp 


. 116 ] 


To determine a and «, values for V 0 , AV, 
Rf, Cf, and w are needed. V 0 and AV are 
measured directly with the measuring lock-in 
amplifier (Fig. 3), a PAR 124. The measuring 
frequency u> is obtained from the front panel 
dial settings of the PAR 1 24, which have been 
checked against calibrated frequency meters. 
The reference frequency w r is measured with 
a digital counter in a Stanford Research Sys- 
tems SR530 lock-in amplifier, which is used 
in the automatic balancing circuit. 

To calculate Rf and Cf, two additional 
measurements or assumptions are required. 
Two quantities are measurable: c T , the con- 
ductivity at the reference frequency, and Cf, 
the reference capacitance. The conductivity of 
the sample at balance is measured using the 
four-electrode cell and a low-frequency con- 
ductivity bridge. Here, the four-electrode cell 
is used in a two-electrode configuration, with 
the source electrodes connected to a bridge 
circuit and the sense electrodes disconnected. 
The cell constant for the chamber used in this 
work is approximately 0.250 cm. An Altex 
(Beckman, Inc.) conductivity bridge is used 
to obtain measurements of the cell conduc- 
tance at 85 Hz and 1 kHz. Since the balance 
frequency is usually in the range 100 Hz to 1 
kHz, these are very convenient frequencies to 
use. Although the conductivity is not mea- 
sured at exactly the balance frequency, mea- 
surements show that the conductivity is not a 
strong function of frequency in this range. 
Thus, the error incurred by the small fre- 
quency mismatch is negligible. 


The reference capacitance, C F , is deter- 
mined as follows. First the capacitors used to 
balance the system are measured with a ca- 
pacitance bridge (General Radio 1620-A) us- 
ing a very accurate three-terminal measure- 
ment at 1 kHz. These capacitances are inde- 
pendent of frequency below 1 MHz. The 
primary error in the measurement of C F is 
that not all of the balance capacitance is di- 
rectly attributable to the measured capacitors. 
For example, stray capacitance to ground at 
the negative sensing lead appears as “extra” 
balance capacitance and is not measurable. 
There are certainly other sources of stray ca- 
pacitance and they cause errors in the deter- 
mination of the effective Cp and the absolute 
dielectric constant, t. The way we deal with 
these errors will be explained shortly. 

Given measured values of o T and Cf, o and 
€ can be computed at any frequency, provided 
that R F can be computed. This is accomplished 
by applying Eq. [ 1 5 ] at the reference frequency 
w r . The output ratio F 0 /AKis written in the 
form 

“7 = r' + />", 117] 

AV 

where r' and r" are both real. At the reference 
frequency they are denoted as r\ and r" so Rf 
is 


Rp = 


1 + r' r 

Go t + o} t C^r'r 


(18] 


Thus, the measurement of oy and Cf allows 
R F to be calculated. Given R f and Cf, a(w) 
and e(w) can be calculated using Eqs. [15] 
and [16]. 

A second approach involves making two 
assumptions instead of two measurements. If 
the conductivity and the dielectric constant of 
the material are known at the reference fre- 
quency, then R f and C F can be computed by 
solving 


1 . . „ G(<j t + /o> r € f Co) 

-R f + '“' Cf " (l+rt + ir?) 


[19] 
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for Ry and C F . This approach is possible with 
pure electrolytes, where a T and c r are known. 

In the first method a T and C F are measured, 
but errors implicit in C F are present. In the 
second approach, the effective value of C F is 
computed exactly, but only for an electrolyte 
where e is known. A hybrid approach can be 
used to take advantage of the favorable attri- 
butes of the two approaches. Before an exper- 
iment with a dilute suspension, measurements 
are made with the pure suspending electrolyte 
and the electrical properties computed by the 
second technique; a T and e T take on their 
“textbook” values at the balance frequency. 
The balance capacitance is then calculated and 
compared to the measured value. A value for 
AC f , the difference between the capacitance 
needed for balance and the measured C F , is 
obtained. Then, in experiments with the sus- 
pension, the first approach is used, but AC F 
is added to the measured balance capacitance 
to compensate for stray capacitance. This pro- 
cedure, which is analogous to the substitution 
method (6) used with conventional bridges, 
reduces the error in the measured values. 

Nevertheless, it should be understood that 
the experiment is not optimized to obtain val- 
ues for a and 6. Instead, it is designed to mea- 
sure the frequency-difference spectra, i.e., <r(v) 
— o(w r ) and c(w) — e(o> r ). The instrument er- 
rors for these quantities are much smaller than 
those for a and e alone. Expressions for the 
difference quantities are 


<r(a>) - <r(o> r ) = 



— U) r i? F C F 


u> 


r - K 


L <o r 


[ 20 ] 


and 


c(a>) - e(w T ) = 


Geo 


1 

<i)r 

— r" - r" 

o) t R f 

a? 


+ C F (r'-r' t ) . 


[ 21 ] 


Consider, first, the expression for a(co) 


- ct(ciy). The most important component on 
the right of Eq. [20] is the first term in the 
brackets, the difference in the in-phase output 
as the frequency is varied between ic T and u>. 
The second term contains the out-of-phase 
output ratio r*. Instead of appearing as a sim- 
ple frequency difference r" — r”, r " is multi- 
plied by the frequency ratio a>/ co r . Since r"(o>) 
is roughly linear in frequency, this term in- 
creases as a; 2 . However, it is multiplied by 
u T RrCr, which is quite small in most circum- 
stances. This prevents the second term from 
being important except at high frequencies. 
Note that only the second term in Eq. [20] 
contains C F so errors in C F do not show up 
until the high-frequency part of the spectrum 
is reached. Tests show that <r(a>) - <r(« r ) at 
high frequencies is indeed sensitive to the value 
of C F . 

The difference spectrum for the dielectric 
constant is much less sensitive to errors in C F . 
According to Eq. [21], the out-of-phase re- 
sponse dominates since the second quantity 
in brackets is multiplied by the small quantity 
C F . The second term contains the frequency 
difference of the in-phase output, which is 
small under most conditions. Thus, errors in 
C F have a small effect on the dielectric con- 
stant difference and experiments show that 
<0) ~ <(o> r ) is insensitive to the exact value 
of C F . On the other hand, the absolute mag- 
nitude of € is sensitive to C F . 

IV. INSTRUMENT CALIBRATION 

In a measurement as complicated as this, it 
is important to provide proof that the results 
are accurate. This is difficult in the range of 
sample conductivities and relaxation frequen- 
cies of interest here since no standards exist 
for a thorough calibration. Nevertheless, our 
testing of the apparatus provides evidence that 
the technique is accurate. 

Two sorts of tests were performed. First, RC 
circuits having a known frequency response 
were substituted for the four-electrode cell. 
The ability of the instrument to duplicate the 
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Fig. 5. The networks used in dummy load tests. 


known response shows that the electronic de- 
sign is sound. Second, “wet” experiments were 
conducted with pure electrolytes of known 
conductivity and dielectric constant, each in- 
dependent of frequency, and with dilute sus- 
pensions. These tests showed that the cell per- 
formed as expected. 

The first electrical tests were designed to 
show that the electronics could reproduce the 
frequency response curve for a“dummy load.” 
Each load was of the general form shown in 
Fig. 5. Terminals “A” and “B”’ represent the 
source electrodes and two 1 00-0 resistors rep- 
resent, in an approximate manner, the imped- 
ance which exists in the regions separating the 
source electrodes from the sense electrodes. 
Terminals “C” and “D” represent the sense 
electrodes. The network between the sense 
electrodes simulates the electrical behavior of 
a material in the cell. 

All of the components used in the test net- 
work were carefully measured: resistances with 
a Data Precision (Analogic Corp.) DVM, and 
capacitances with a capacitance bridge at 1 
kHz. Tests were carried out with five different 
loads to cover the range of responses expected 


with the cell. In each test an experiment was 
conducted with a given dummy load in place, 
and the in-phase and out-of-phase compo- 
nents of the ratio 1 o/AFwere measured. Then 
the predicted frequency response was com- 
pared with the experimental results. Figure 6 
shows results from one such test. Agreement 
between the prediction and the experiment 
was always excellent. 

Tests were also carried out with aqueous 
electrolytes. Since the conductivity and di- 
electric constants of these materials are inde- 
pendent of frequency over the frequency range 
of interest, the measured frequency-difference 
spectrum should be flat. Electrical tests showed 
that despite the good agreement between 
measured and predicted output voltages, there 
is an electronic limitation to the accuracy of 
the dielectric constant measurement at low 
frequencies. Electrolyte tests show whether 
electrochemical problems cause additional er- 
rors. 

Results from one set of tests are shown in 
Fig. 7. In the upper panel the absolute con- 
ductivity is plotted against frequency. Note 
that the scale is expanded considerably; the 
conductivity is essentially constant to 20 kHz, 
and the error at 200 kHz is less than 2%. The 
second panel shows the dielectric constant. 
The percentage error shown for this test is 
about 20%. In most cases the error in the di- 
electric constant is much smaller, as will be 
shown shortly. Furthermore, the error de- 
creases as the dielectric constant of the ma- 




Fig . 6. Comparison between measured and predicted output ratios for a dummy load test showing the 
in-phase (left) and out-of-phase (right) data; theory, — ; experiment, O. 
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Fig. 7. Measured conductivities and dielectric constants 
for a 10~ 4 M HC1 solution. 


terial in the cell increases so that the errors in 
the dielectric behavior of suspensions are 
much less [6]. 

V. TESTS WITH SUSPENSIONS 

Extensive tests were carried out to establish 
the characteristics of the instrument with sus- 
pensions in the measuring cell. The suspen- 
sions consisted of 0. 1 9 1 -jtm amphoteric latex 
particles in 10 -4 M HQ. Detailed information 
on these particles is given in a companion pa- 
per ( 9 ) . For the present it suffices to note that 


the ^-potential of these particles in 10" 4 M 
HCl is roughly 80 mV. To obtain frequency- 
difference spectra, six separate experiments 
were conducted, two with the pure electrolyte 
and four with suspensions. The volume frac- 
tion 0 varied between 0.024 and 0.04 1 in the 
experiments with suspensions. The reference 
frequency used in every case was 500 Hz 
and the measuring frequencies, eight in all, 
followed a logarithmic sequence from 1 to 
200 kHz. 

Figure 8 shows the variation of the absolute 
dielectric constant with frequency at different 
particle concentrations. Results of the baseline 
measurements before and after the tests with 
suspension are almost identical, and flat. The 
repeatability, and the minimal frequency de- 
pendence of the measured dielectric constant 
of the pure electrolyte, indicates that no un- 
usual electrochemical effects influenced the 
experiment. Note also that the accuracy here 
is better than that found in the earlier tests. 

The remaining curves show results for ex- 
periments with suspensions. As expected, there 
is a tremendous change in the dielectric con- 
stant due to the addition of relatively small 
amounts of colloidal material; a dielectric 
constant of 800 with a volume fraction of only 
0.041 is 10 times the value for pure water. 
Suspensions with successively smaller particle 
volume fractions yield smaller values for the 
low-frequency dielectric constant. In each case 
the dielectric constant decreases as the fre- 
quency is increased. At 50 kHz, the dielectric 
constant of the suspension is very close to that 
of pure water. 

To further check the instrument, the fre- 
quency difference spectra measured with sus- 
pensions of various particle concentrations 
were normalized by the volume fraction of 
particles and plotted together (Figs. 9 and 10). 
For dilute suspensions the spectra should be 
linear in the volume fraction and the fact that 
the data superimpose shows that the cell and 
electronics are free of errors due to the elec- 
trical or chemical properties of the suspension. 
If such errors did exist, they would need to be 
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FlG. 8. The dielectric constant for suspensions of (0.191-^m) latex parties in 10' 4 M HQ; volume 
fractions: O, 0.041; □, 0.035; A, 0.032; 0 , 0.024; ®, 0.000. 



Fig. 9. Normalized frequency-difference spectra for the dielectric constant of suspensions of (0.191-mn) 
latex particles in 10" 4 M HC1; volume fractions as noted on Fig. 8. 
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Fig . !0. Normalized frequency-difference spectra for the conductivity of suspensions of (0. 191-^m ) latex 
particles in 1CT 4 M HC1; volume fractions as noted on Fig. 8. 


linear in the volume fraction of particles to 
avoid detection, which is unlikely. 

The linearity of the frequency-difference 
spectra with respect to particle volume fraction 
was also checked at every frequency using lin- 
ear regression. The linear correlation coeffi- 
cient for the experiments was in excess of 0.998 
except for the conductivity measurements at 
100 and 200 kHz. 

The change in the conductivity caused by 
the addition of particles to the electrolyte is 
markedly less than the change in the dielectric 
constant. The conductivity of 10 ~ A M HC1 is 
about 45 ftS/cm, and a 4.1% by volume ad- 
dition of particles increases this by about 2 
/xS/cm. Instead of an order of magnitude 
change, this is a modest 4.4% increase. This 
makes the conductivity frequency-difference 
spectrum more difficult to measure than the 
dielectric constant difference spectrum. 


The experiments described were repeated 
using a fresh suspension as a check on the re- 
producibility of the results. The differences in 
the frequency-difference quantities were less 
than 5%. 

VI. DISCUSSION 

The instrument described here employs a 
four-electrode frequency-difference technique 
adapted for colloidal suspensions. Dielectric 
spectrometers of this type are an improvement 
over the traditional Wheatstone bridge devices, 
offering shorter measurement times, improved 
signal-to-noise ratios, and elimination of the 
thermal drift problem without precise tem- 
perature control. The dielectric constants and 
conductivities of colloidal dispersions mea- 
sured with this instrument are linear in the 
volume fraction of particles at low volume 
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fractions and quite reproducible. In the sub- 
sequent paper (9) we discuss extensive mea- 
surements covering a range of ionic strengths. 
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A four-electrode frequency-difference dielectric spectrometer is used to study the dielectric response 
of suspensions of amphoteric polystyrene particles in dilute HC1 solutions. Results are given for one set 
of particles with relatively thick diffuse layers, i.e., 3 < ax < 10. When the dielectric spectra are compared 
to those predicted by classical electrokinetic theory, large discrepancies are found. The measured dielectric 
constant is an order of magnitude larger than the theoretical value. The limited body of experimental 
data extant in the literature confirms these results, suggesting that the standard electrokinetic model is 
somehow deficient when applied to dispersions of this sort. © 1989 Academic Press, Inc. 


I. INTRODUCTION 

A growing body of evidence suggests that 
the electrokinetic behavior of some colloidal 
particles cannot be characterized by a single 
parameter such as the ^-potential (1-4). 
However, much of the experimental evidence 
is based on static measurements, such as the 
electrophoretic mobility of individual parti- 
cles. or low-frequency measurements (<100 
Hz ) , such as the electrical conductivity of sus- 
pensions. Thus, to some extent, the evidence 
is inconclusive. Part of the difficulty stems 
from the low information content of these 
measurements, essentially one data point per 
experiment. Furthermore, since the anoma- 
lous behavior appears connected with pro- 
cesses on the particle surface, it is important 
to be able to study this region and its inter- 
action with the diffuse layer experimentally. 
Dielectric spectroscopy offers a means of 
studying relaxation behavior because the 
properties of a suspension can be measured 
over a range of frequencies. In the first paper 
of this series a new instrument was described 


1 To whom correspondence should be addressed. 


(5). Here we report measurements on disper- 
sions made at several salt strengths. 

II. EXPERIMENTAL PROTOCOLS 

Suspensions of polystyrene particles are 
popular model colloidal systems because 
standard synthesis techniques produce nearly 
spherical and monodisperse particles. An ad- 
ditional advantage is that surface properties 
and the particle diameter can be varied by 
modifying the synthesis procedure. An am- 
photeric polystyrene latex served as the model 
system for this work. It was synthesized ac- 
cording to the procedure developed by Hom- 
ola and James (6 ), with the modification sug- 
gested by Harding and Healy ( 7 ) . Particle sur- 
faces contained both amine and carboxyl 
groups. The former are positively charged at 
low pH and neutral under other conditions; 
the latter are negatively charged at high pH. 
On the basis of the synthesis used for our latex, 
the isoelectric pH should be approximately 5.4 
(7). In our experiments the particles were in 
the cationic form. 

Particle size and monodispersity were de- 
termined by transmission electron micro- 
scopy, which gave an average particle diameter 
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of 0.191 with a standard deviation less 
than 0.003 jum. Particles in the TEM images 
appeared to be perfectly spherical. 

The particles were cleaned using a sedi- 
mentation-decantation procedure (4) that 
employed a Beckman Model L5 ultracentri- 
fuge with the Model 21 rotor running at 8000 
rpm (9000 g). A run of 90 min was sufficient 
to achieve complete separation and produced 
a pellet which was easily redispersed without 
sonication. Particles were spun down in elec- 
trolyte at the ionic strength used in the exper- 
iment. The cleaning process continued over 
several sedimentation-decantation cycles until 
the supernatant conductivity was within 1% 
of the conductivity of the pure suspending 
electrolyte. The spectroscopy experiment was 
completed within 48 h of cleaning. 

Particle mobilities were measured by mi- 
croelectrophoresis, using a Rank Brothers 
Mark II device. Potentials were computed 
from the measured electrophoretic mobilities 
with the computer program of O’Brien and 
White (8); results are shown in Fig. 1. These 
data represent two series of experiments, con- 
ducted approximately 30 days apart with the 
error bars reflecting two standard deviations 
about the mean of 20 measurements at each 
salt strength. The large error bars near 0.001 
M reflect the flat nature of the f-potential- 



Fig. 1. The f-potentials of 0.191-pm particles in HC1 
solutions. 


mobility relation at this double-layer thickness; 
the scatter in the mobility measurement is 
much less. 

Prior to each experiment involving a sus- 
pension, the cell used in the dielectric spec- 
troscopy experiment was filled with the sus- 
pending electrolyte and allowed to stand 
overnight. Then, immediately before the ex- 
periment, the cell was emptied and refilled 
with electrolyte approximately five times. A 
sequence of six experiments was performed 
for each ionic strength. First, an experiment 
with the pure suspending electrolyte was done 
to establish the experimental baseline. Next, 
the cell was filled with suspension and four 
experiments were conducted, each at a differ- 
ent particle volume fraction. Finally, the cell 
was emptied of suspension and thoroughly 
cleaned with the pure electrolyte, and the 
baseline measurement was repeated. 

The experiment sequence proceeded toward 
successively less concentrated suspensions, 
each obtained by withdrawing a small volume 
of suspension with a pipet and replacing it with 
the same volume of the pure electrolyte. Then, 
the contents of the cell were mixed by redis- 
persing with a pipet. The small volume ( < 1 
cc) of suspension removed during each dilu- 
tion was saved for the volume fraction deter- 
mination. Each sample was saved in a plastic 
vial and refrigerated until the entire set of ex- 
periments was complete. The volume fraction 
was measured by drying a known mass of sus- 
pension. assuming the specific gravity of poly- 
styrene to be 1.057. 

Several sets of experiments were done with 
particles suspended in HC1 at increasing salt 
strengths. The particles for the first experiment 
( 1G -4 M HC1 ) were cleaned starting from raw 
latex. Subsequent experiments used the same 
latex cleaned in the more concentrated salt 
solution, with a small quantity of the raw latex 
stock added to make up for the sample used 
in the volume fraction determination. 

III. EXPERIMENTAL RESULTS 

The experimental results were analyzed us- 
ing procedures described by Myers and Saville 
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( 5 ) . Experiments on suspensions with different 
particle concentrations were checked to ensure 
that the response was linear in the volume 
fraction of particles. Only data for which the 
correlation coefficient was better than 0.99 are 
reported; in most cases the correlation coef- 
ficient exceeded 0.999. The results are reported 
as “difference spectra,” i.e., the difference be- 
tween the value at the measuring frequency 
and the reference frequency (500 Hz) nor- 
malized by the volume fraction of particles, 
<t>, i.e., 

conductivity spectrum: [a(a>) — <r r ]/<£ 

dielectric constant spectrum: [e(a>) - e r ]/<£. 

Note that this procedure eliminates effects due 
to nonspecific adsorption and added counter- 
ions ( 9 ) . 

Figure 2 shows the normalized frequency- 
difference spectrum for the conductivity at 
several ionic strengths. The most striking fea- 
ture is the increase in the magnitude of the 


conductivity difference with increasing salt 
strength. Increasing the HC1 concentration by 
a factor of 10 brings a factor of 4 increase in 
the conductivity difference at 50 kHz. The full 
behavior of the relaxation frequency is not 
readily apparent here. The data at 10 4 M 
reach the high-frequency plateau at 50 kHz, 
but at the other ionic strengths there is still a 
substantial frequency dependence at the high- 
est frequencies measured. 

The normalized frequency-difference spec- 
tra for the dielectric constant are shown in Fig. 
3. Here again, the magnitude of each spectrum 
increases with increasing salt strength. The di- 
electric constant difference has clearly reached 
'a high-frequency plateau by 20 kHz in 10~ 4 
M HC1, but not with 1 0 ~ 3 M HC1. At inter- 
mediate salt strengths the same trend is ap- 
parent, leading to the conclusion that the 
characteristic frequency increases monotoni- 
cally with salt strength. 

The experiments described in this section 
were repeated with fresh latex as a check on 



Fig. 2. Frequency difference spectra for the conductivity of suspensions of 0. 191-jtm particles at different 
salt strengths; O, 10~ 4 M, uk — 3.1; D , 2 X 10 4 M, ar — 4.4; A, 5 X 10 A/, cik — 7.0, v , 10 A/, cik 

= 9.9. Note that the frequencies quoted here and elsewhere in the figures are the linear frequencies. w/27t. 

Journal of Colloid and Interface Science. Vol. 131, No. 2. September 1989 



464 


MYERS AND SAVILLE 



Fig. 3. Frequency difference spectra for the dielectric constant of suspensions of 0. 191-^m particles at 
different salt strengths; O, 10 -4 M, ok = 3.1; □, 2 X 10 -4 A/, an = 4.4; A. 5 X 10 -4 A/, an = 7.0; v, 10 -3 
M, an = 9.9. 


reproducibility. The average deviation for the 
conductivity was 7%; for the dielectric constant 
it was 4.6%. 

These data show that the conductivity and 
dielectric constant frequency-difference spec- 
tra are strong functions of ionic strength, in 
terms of both the magnitude of the effect and 
the relaxation frequency of the dielectric con- 
stant. The spectra are linear in the particle 
volume fraction and reproducible. In the next 
section these results are compared with a the- 
ory of dielectric relaxation based on the clas- 
sical electrokinetic model. 

IV. A COMPARISON BETWEEN 
EXPERIMENTAL DATA AND THEORY 

Since bulk suspension properties alone are 
accessible to our experimental technique, a 
mathematical model must ‘ be used to relate 
the relevant physical processes to the surface 
properties of the particles. Here, the experi- 
mental data are compared to predictions from 
the standard electrokinetic model, using the 
solution of DeLacey and White (10). Al- 
though a number of dielectric relaxation 


models for colloidal suspensions have been 
proposed in the literature (11-16), only the 
work of DeLacey and White (10) treats the 
problem without any mathematical approxi- 
mations. Theirs is a thorough treatment of the 
problem using the standard model. To avoid 
ambiguities, the equations making up the 
standard model will be reviewed briefly. 

In the classical theory of electrokinetics, 
particles of radius a and dielectric constant e p 
are immersed in an electrolyte containing N 
distinct ion types. Each ion is characterized 
by its volume density rij, mobility m } , and va- 
lence Zj. The electrolyte is assumed to act as 
a Newtonian fluid with viscosity p, density p, 
and dielectric constant . N + 3 coupled dif- 
ferential equations and N + 1 constitutive re- 
lations are available to solve for the indepen- 
dent variables: the fluid velocity u; the fluid 
pressure P\ the electric potential \p; and the 
ion densities j = 1,2,..., N. For oscil- 
latory fields, the time-dependent nature is as- 
sumed to be exp(zW); ^ is the angular fre- 
quency of the applied field and z stands for 
V— T. The differential equations are 
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€o«oo vfy = -p e [1] 

/^V 2 u — VP = — p e E [2] 

V • u = 0 [3] 

j=l,2,...,N [4] 

N 

Pe = Z l 5 ] 

V=1 

fj - riju - m } ez } np^ — mjkTVnj. [ 6 ] 


Equation [1] is Poisson’s equation, with p e 
representing the volume charge density in the 
electrolyte. Equations [2] and [3] are the 
Stokes equations where viscous and pressure 
forces are balanced against the electrical body 
force term; E is the electric field vector. Equa- 
tion [4] relates the divergence of the flux of 
the yth ion to the time rate of change of its 
volume density. The constitutive relation for 
charge density is given by Eq. [5], where e is 
the fundamental electric charge. Equation [6] 
gives the ion fluxes in terms of components 
due to convection, electromigration, and dif- 
fusion. Here k is Boltzmann’s constant and T 
is absolute temperature. The boundary con- 
ditions will be discussed shortly. 

It is desirable to separate this problem into 
two simpler problems by assuming that the 
dependent variables can be written as a per- 
turbation expansion in the dimensionless 
electric field strength (3 = aeE 0 /kT, where E 0 
represents the imposed field. Thus, for (3 1, 

yp = + W + 0((3 2 ) 

rij = n® + (38rij + 0((3 2 ). [7] 

First, the 0{ 1 ) problem, which represents 
the equilibrium situation, is solved with the 
Poisson-Boltzmann equation used to relate \j/° 
to rij . The boundary conditions on the electric 
potential are 

\p°(r) —*• 0 as | r | ->■ cc [8] 

PVW- [9] 

The 0(/3) problem is more complicated. 
Fluid flow must be considered and, due to the 


oscillatory nature of the problem, all the terms 
in the perturbation equations are complex- 
valued. Boundary conditions must be applied 
to 8\p, <5u, and all of the 8nj in the O(fi) prob- 
lem. Far from the particle, 


d\p -*• — E 0 • r 

[10] 

5u -*■ — /i e E 0 

[113 

brij -*■ 0, 

[12] 

where ii e stands for the electrophoretic 
bility of the particle. 

At the particle surface, 

mo- 

fy(a) - 8Ma) = 0 

113] 

[tV8\j/(a) — e p V5i p\(a)] • n = 0 

[14] 

n = 0 

[15] 

u(tf) = 0. 

[16] 


Equation [13] ensures that the electric poten- 
tial is continuous across the particle surface 
(8\pi is the electric potential inside the particle). 
The next expression represents the charge in- 
duced at the particle, which is proportional to 
the jump in the electrical displacement; n is a 
unit vector pointing out from the surface. 
Since the jump is set equal to zero, the particles 
have a constant charge. The other equations 
specify that the ion fluxes normal to the par- 
ticle surface vanish and that the no-slip con- 
dition prevails. 

These model equations constitute what is 
called the standard model of electrokinetics: 
the Gouy-Chapman model of the diffuse layer 
coupled with conservation equations and 
boundary conditions where the equilibrium 
interface is unaffected by dynamic processes. 
DeLacey and White (10) solved these equa- 
tions numerically to obtain the dipole coeffi- 
cient from the far-field form of the electric po- 
tential and this was used to compute the com- 
plex conductivity, i.e., the conductivity and 
dielectric constant of the suspension. Their re- 
sults are presented in a normalized form which 
is independent of the volume fraction of par- 
ticles using two functions; Ae' and A/f efr . These 
are defined by 
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e = toe + 0Ae' [17] 

<r = tr,* + 0A^ eff , [18] 

where and are the dielectric constant 
and conductivity of the electrolyte in the ab- 
sence of particles, and 0 is the particle volume 
fraction. 

The frequency-difference spectra from our 
dielectric spectroscopy measurement can be 
related to DeLacey and White’s (10) functions 
in a straightforward way using Eqs. [17] and 
[18]. Thus 

<r(aj) -- - = AK eff (co) - AK eff (u> r ) [19] 

0 

€(C0) ~ - = A<*w) - Ae(w r ). [20] 

0 

Accordingly, the normalized frequency-dif- 
ference spectrum from the measurement is 
equal to the frequency-difference spectrum of 
DeLacey and White’s functions. 

DeLacey and White’s published results (10) 
have been checked for low ("-potentials by 
O’Brien (17), who found good agreement with 
dimensionless ("-potentials as high as 2. The 
DeLacey-White algorithm is available as a 
FORTRAN program, which we have imple- 
mented and checked by computing the results 
shown in three figures in DeLacey and White’s 
paper (Figs. 1, 4, and 7, each at f = 100 mV). 
We also checked the zero-frequency limits of 
the conductivity and mobility against the re- 
sults of Zukoski and Saville (18). 

The first comparison with experimental 
data will be for dispersions in 10“ 4 M HC1. 
The measured ("-potentials are given in Fig. 1. 
Figure 4 shows the measured and theoretical 
frequency-difference spectra for the conduc- 
tivity. The theoretical spectrum is about 50% 
larger than the measured magnitude but the 
characteristic relaxation frequencies differ by 
over an order of magnitude. 

The comparison between measured and 
theoretical spectra for the dielectric constant 
shows even larger discrepancies ( Fig. 5 ) . Here 
the predicted magnitude of the relaxation is 
more than an order of magnitude smaller than 



Fig. 4. Comparison between the measured frequency- 
' difference spectrum for the conductivity and that predicted 
by the standard model for 0.19 1-^ni particles in 10 ~ 4 M 
HC1. 

that measured. The difference between the re- 
laxation frequencies is similar to that observed 
in the conductivity-difference spectrum. Sim- 
ilar differences were found at the other salt 
concentrations ( Figs. 6 and 7 ) . 

Needless to say, these are very striking 
comparisons. The discrepancy between theory 
and experiment seen in the conductivity spec- 
trum magnitude is relatively small and might 
be brought into agreement by adjusting the ("- 
potential, but the predicted relaxation fre- 
quency is far too large. However, both the 
magnitude and the characteristic frequency of 
the theoretical dielectric constant spectrum 
disagree considerably with the experimental 
results. These discrepancies cannot be re- 
moved by adjusting ("-potentials. 

~ "Several factors were investigated to see if 
they would account for the discrepancies: er- 
rors in the ("-potential and the formation of 
permanent or temporary aggregates. In the 
experiments at 10' 4 M HC1 the measured ("- 
potential was 100 mV. To test the sensitivity 
of the theory to changes in the ("-potential, 
results were calculated for ("-potentials of 100 
and 200 mV. Figure 8 shows that this effect is 
too small to account for the discrepancy. 

Another possibility is the presence of ag- 
glomerates. The predicted dielectric constant 
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FREQUENCY (KHz) 


FiG. 5. Comparison between the measured frequency- 
difference spectrum for the dielectric constant and that 
predicted by the standard model for 0. 19 l-jun particles in 
10" 4 A/HCI. 


FiG . 7. Comparison between the measured frequency- 
difference spectrum for the dielectric constant and that 
^predicted by the standard model for 0. 1 9 1 -/im particles at 
four salt strengths: O, 10 -4 M, cu = 3.1: □. 2 X 10" 4 M, 
ok = 4.4: A, 5 X 1(T 4 Af, ok = 7.0; 0, 10" 3 M. ok = 9.9. 


increases with particle size and Lyklema et ah 
(19) suggested that this could lead to discrep- 
ancies. Increasing the value of an by as much 
as a factor of 16, however, fails to bring theory 
and experiment into agreement and it is un- 
likely that agglomerates of this size or larger 
existed in our suspensions. No agglomerates 
of any size were seen in the TEM images, so 
the agglomeration would have to take place 
in the cleaning procedure or during the ex- 


periment. Furthermore, agglomerates of this 
size would have settling properties different 
from those of individual particles. No sedi- 
ment was ever observed in the bottom of the 
sample cell, nor was any coagulum ever ob- 
served floating on the surface of the sample. 
It can be stated with some confidence that ag- 
glomerates of this size did not exist before or 
after our experiments. 



Fig. 6. Comparison between the measured frequency- 
difference spectrum for the conductivity and that predicted 
by the standard model for 0. 19 l-£tm particles at four salt 
strengths: O. 10~ 4 M, ok = 3.1; □. 2 X 10 4 M , ok = 4.4; 
A, 5 X 10' 4 A/, ok - 7.0; >, 10“ 3 M, ok = 9.9. 



Fig. 8. Comparison between the measured frequency- 
difference spectrum for the dielectric constant and that 
predicted by the standard model for 0. 1 9 1-Atm particles in 
10~ 4 M HC1 at two f-potentials. 
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A third possibility is that temporary aggre- 
gates were formed. If dipole-dipole interac- 
tions between particles induced doublets or 
larger chains to form, this would increase the 
dielectric constant and the doublets would 
disappear after the field was removed. How- 
ever, this mechanism can be discounted on 
the grounds of a scaling argument and exper- 
imental evidence. The interaction energy be- 
tween two dipoles separated by a distance r is 
given by Israelachvili ( 20 ) as 

2u\U 2 

Interaction Energy = t , [21 ] 

47 t te 0 r 

where U\ and u 2 are the dipole moments of 
the interacting particles. If doublets are 
formed, the interaction energy must exceed 
the energy of thermal motion, kT. Thus, for 
identical particles dipole-dipole interactions 
will be negligible when 


2 u 2 

AirttQr^kT 


< 1. 


[ 22 ] 


To compute this ratio, an estimate for the di- 
pole moment is needed. This is related to the 
particle dipole coefficient, i.e., 

u = 4xeeo£ 3 £oCo> [23] 


with C 0 representing the dipole coefficient of 
a particle. Thus, 


64xec 0 Q)-^o^ JI 

kT 


<§ 1 


[24] 


if dipole-dipole interactions are negligible. 
Here, the interaction distance r has been as- 
sumed to equal to twice the particle radius a. 
Taking C 0 = 1.0 (a large value), an electric 
field strength of 1 0 V / m ( typical of those used 
in our experiments), and a particle radius of 
0.1 jum, yields a value of 3.6 X 10~ 6 . This is 
sufficiently small to justify neglecting dipole 
interactions. 

Finally, to see if the field strength in our 
apparatus had any effect, frequency scans were 
done with field strengths of roughly 10, 20, 
and 40 V/m. Since the scaling argument in- 
dicates that the dipole-dipole interaction in- 


creases with the square of the electric field 
strength, doubling the field should increase the 
dipole interaction by factors of 4 and 16. Little 
or no effect was observed and we conclude 
that effects due to induced dipoles are negli- 
gible. 

At this point we are forced to conclude that 
the standard model fails to describe our ex- 
perimental results. Next, other experimental 
studies were scrutinized to see whether simi- 
larly large deviations between experiment and 
theory had been observed previously. The 
work of Springer et al. (21) and Lim and 
Franses (22) was used although precise ("-po- 
tential data were not available. Springer et al. 
do not report ("■ -potentials, but do give titrat- 
able charge from which we have estimated the 
('-potentials. Lim and Franses report ("-poten- 
tials, but they have an unexplained depen- 
dence on the particle volume fraction and may 
be inaccurate. Although the ("-potentials may 
be in error, the influence of the errors can be 
estimated. It should also be noted that neither 
of these reports included a direct comparison 
of experimental data with the standard model 
in a detailed form. However, Springer el al. 
did compare their results with a theory for thin 
double layers and found substantial disagree- 
ment (19). 

For the purposes of comparison, a measure 
of the magnitude of the change in the dielectric 
constant is needed. We used the “static specific 
dielectric increment” (SSDI) defined by Lim 
and Franses (22), viz., 


SSDI = 


*(0) - e(oo) 


[25] 


where 

«(0) = lim e(cij) and 


e(oo) = lim do;). [26] 

a>— ►ac 

In each case studied, the dielectric constants 
in the two limits and the particle volume frac- 
tion were obtained from the published exper- 
imental data and computed with DeLacey and 
White’s computer program (10). For com- 
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parison with Lim and Franses’ work the mea- 
sured ("• -potentials were used; for Springer et 
al. (21 ), ("-potentials were estimated. 

It should be pointed out that there is some 
ambiguity in the choice of t(0) and e(oo) in 
each data set. For Lim and Franses’ data (22), 
the published values based on a Cole-Cole ex- 
trapolation were used. For Springer et al. (2 1 ), 
as well as our own data, the lowest frequency 
and highest frequency data points were used. 
Although these procedures are imprecise, we 
are looking for order-of-magnitude effects so 
the approximation is useful. The SSDIs from 
the experiments will be underestimated, in any 
event. 

The results of the comparison are given in 
Table I. To test the sensitivity to changes in 
the ("-potential, different values were chosen. 
Nevertheless, the differences cannot be ex- 
plained by low ("-potentials. The tabulated data 
show that at least three different experimental 
studies yield dielectric constants substantially 
larger than predicted by the standard model. 
Moreover, a large number of experimental 
parameters have been varied in the studies. 
Both Springer et al. (21) and Lim and Franses 
(22) used two-electrode bridge methods, while 
a four-electrode frequency-differencing tech- 


nique was used in our work. It seems unlikely 
that experimental errors would consistently 
give inflated values for the dielectric constant 
using such different techniques. 

V. CONCLUSIONS 

A comparison between the experimental 
results and the predictions of the standard 
model has identified several discrepancies. The 
difference between theory and experiment in 
the measurement of the electrical conductivity 
is relatively small, consistent with the earlier 
results of Van der Put and Bijsterbosch (1,2), 
O’Brien and Perrins (3), and Zukoski and Sa- 
ville (4). However, the characteristic relaxa- 
tion frequency of the suspension and the mag- 
nitude of the dielectric constant, two mea- 
surements unique to dielectric spectroscopy, 
show substantial differences between theory 
and experiment. Studies with other instru- 
ments show similar behavior. 

Our work and the work of other investiga- 
tors involve a variety of configurations and 
conditions: two- and four-electrode cells, two 
types of particles ( polyvinvltoluene and poly- 
styrene), three different salts ( NaCI. KCI, and 
HC1), and a range of salt strengths and particle 


TABLE I 


Comparison between Theoretical and Measured SSDI Values 


Case 

Author 0 

Salt 

a (jim) 

an 

f (mV) 

SSDI X I0' z 
(experimental) 

SSDI X I0' J 
(theoretical) 

1 

L&F 

NaCI 

1.0 

120 

-37 

1300 

4 (307)* 

2 

L&F 

NaCI 

0.6 

80 

-33 

2000 

3 (630) 

3" 

L&F 

NaCI 

0.6 

80 

-132 

2000 

250 (8) 

4 

S 

KCI 

0.2 

15 

-50* 

170 

7 (23) 

5 

S 

KCI 

0.2 

15 

-75* 

170 

20 (9) 

6 

S 

HC1 

0.2 

12 

1 

o 

140 

8(18) 

7 

S 

HC1 

0.2 

12 

o 

o 

T 

140 

41 (3) 

8 

M&S 

HC1 

0.1 

3 

81 

170 

15 (11) 

9 

M&S 

HC1 

0.1 

4 

84 

200 

18 (11) 

10 

M&S 

HC1 

0.1 

7 

103 

250 

32 (8) 


“ L&F, Lim and Franses (22): S. Springer el al. (21); M&S. this work. 

* Numbers in parentheses are the ratios of the experimental to theoretical SSDIs. 
r Conditions the same as for case 2 except that the ('-potential has been quadrupled. 
d Estimated ("-potentials for each case. 
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sizes (an varied from 3 to 120). In every case 
the experimental data on the dielectric con- 
stant exceeded the theoretical values obtained 
from the standard model. This suggests that 
the discrepancies encountered with synthetic 
latex particles result from a general phenom- 
enon somehow absent from the standard 
model. 
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ABSTRACT 


Recent developments in theory and instrumentation have contributed to renewed interest in the 
dielectric properties of colloidal dispersions. One reason for the renewed interest is that dielectric 
measurements appear to be more sensitive to characteristics of the particle surface than conventional 
electrophoretic mobility measurements. However, the amount of data on well-characterized systems is 
small and it is premature to generalize. To broaden the data base, two different types of colloidal 
particles were studied: an amphoteric polystyrene latex and colloidal silica. In both cases, the 
dielectric response was found to differ from that predicted with the standard electrokinetic theory using 
^-potentials measured electrophoretically. The dielectric constant of the latex suspensions exceeded 
the theoretical prediction by an order of magnitude or more. At the same time, the characteristic 
relaxation frequencies were consistently lower than the theoretical. Experiments with silica 
dispersions yielded smaller dielectric constants than those found with the latex and the disagreements 
between theory and experiment were less. The results also indicate that dielectric spectroscopy is 
sensitive to characteristics of the dispersion other than the particle size and the ionic composition of 
the suspending electrolyte. 
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INTRODUCTION 


The electrostatic state of colloidal particles is often characterized by a single quantity, the £- 
potential. However, for polystyrene latiees this parameterization can be ambiguous since ^-potentials 
measured by different techniques often disagree 1 ' 3 . This is indicative of a fundamental problem since the 
classical electrokinetic theory 4 ' 6 is used to interpret the measurement in each study. It has been 
suggested that the inconsistency in the ^-potential results from disregarding processes occurring on the 
particle surface 7 . Before this hypothesis can be tested comprehensively, techniques need to be developed 
which yield more information per experiment than mobility or static conductivity measurements. 
Dielectric spectroscopy offers an intriguing possibility, but our experience with well-characterized 
suspensions is small. The purpose of this work is to increase our understanding of the dielectric behavior 
of model dispersions. 

Dielectric spectroscopy yields much more information than static measurements of mobility or 
conductivity due to the presence of dielectric relaxation at frequencies below 100 kHz 8 ’ 9 . This relaxation 
is of particular interest since it arises from the temporal behavior of transport processes which polarize 
the electrical double layer 10 . The effects due to diffuse layer transport are included in the classical 
electrokinetic theory 4 ’ 11 , and predictions for the dielectric constant and the conductivity are available 
as a function of frequency. Thus, given the particle size, the ^-potential, and the electrolyte 
composition, experimental results can be compared with the theory. 

Previous measurements of dielectric spectra for anionic latiees 9,12 and an amphoteric latex in the 
cationic form 13 reveal that the low-frequency dielectric constant is substantially larger than the 
classical theory predicts. Furthermore, increasing the ^-potential to values much greater than those 
inferred from the electrophoretic mobility does little to improve matters 13 . Evidently, the small, 
systematic discrepancies in the static-field experiments 1 ' 3 are greatly amplified in measurements of the 
low-frequency electrical response. In this paper, we present results for positively-charged polystyrene 
particles and negatively-charged silica spheres. The objectives are to examine the effect of the 
counterion on the suspension's electrical response, and to expand the experimental database in 
anticipation of future modelling work. Measurements of the normalized frequency-difference spectra are 
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reported for latex in two different electrolytes (HC1, HNO3) with 4 < aK < 12. Electrophoretic 
mobilities were measured to facilitate a comparison of the dielectric data with the classical 
electrokinetic theory. It is found, once again, that the magnitude of the relaxation is much larger than 
expected by the theory, and that the characteristic relaxation frequency is lower than the theory 
predicts. The dielectric response for the silica systems differs markedly from that of the latex, 
suggesting that the electrical properties of colloidal suspensions are sensitive to the specific nature of the 
particle surface. 


EXPERIMENTAL 

Amphoteric Latex 

An amphoteric polystyrene latex was synthesized according to the procedure of Homola and 
James 14 , using the modification discussed by Harding and Healy 1 ^. The amphoteric (or zwitterionic) 
latex possesses amidine and carboxyl groups which give rise to a pH-dependent surface charge. Based on 
the synthesis recipe, the predicted isoelectric pH of the latex is approximately 5.4. For our experiments, 
the latex was in the cationic form. 

The latex was cleaned by sedimentation/decantation in a Beckman Model L5-65 ultracentrifuge 
with a Type 21 rotor. Dilute hydrochloric acid was used initially as the suspension medium, followed by 
the electrolyte to be employed in the electrokinetic experiments. Approximately twelve cycles were 
required to remove residual ionic impurities. The cleanliness of the latex was assured by requiring that 
the decanted supernatant conductivity be constant to within 1% of the conductivity- of the wash 
electrolyte. Experiments were carried out in sequences of increasing ionic strength, beginning with O.lmM 
HC1 and culminating with ImM HC1. Latex equilibrated in O.lmM HQ was used as the starting point for 
the nitric acid experiments, beginning with O.lmM HNO3. 

Particle size and monodispersity were determined by transmission electron microscopy. The 
magnification was calibrated with latex standards of known diameter (Interfacial Dynamics, Portland 
OR). The latex had a mean diameter of 226 ran with a coefficient of variation of 3.0%. 
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Electrophoretic mobilities were measured with a PenKem 3000 Instrument. The HC1 samples 
inexplicably showed two peaks corresponding to small positive and negative mobilities (< 1 pm-cm/volt- 
sec in absolute magnitude), suggesting that those samples were somehow contaminated. The HNO3 
samples had uniformly large positive mobilities (> 4 |im-cm/volt-sec) which were comparable to those 
measured previously for a similar latex 13 . Accordingly, we used the HNO3 mobility data in the 
comparisons with theory. Mobilities for the HNO3 samples (Table I) were converted into ^-potentials 
using the computer program of O'Brien and White 4 . 

Silica 

Colloidal silica furnished by E. Matijevic was prepared according to the method of Stober et al. 16 
with a proprietary modification used to synthesize large quantities of material 17 . The particles were 
aged in dilute electrolyte for approximately two months, whereupon the pH was adjusted to 4.0 
with 0.1M HC1. The dispersion was equilibrated at the desired pH and ionic strength by 
sedimentation/decantation, using six wash cycles over roughly a two week period until the conductivity 
and pH of the supernatant were constant. Once again, dielectric spectroscopy measurements were made in 
an order of increasing ionic strength. 

The mean particle diameter was 330 nm with a coefficient of variation of 7.2% according to 
transmission electron microscopy. Zero-field electrophoretic light scattering with the Coulter DELSA 
440 instrument provided a rough estimate of the mean diameter, 300 nm < d < 350 nm. 

Electrophoretic mobilities were measured with the Coulter DELSA 440 (Table I). ^-potentials were 
obtained from the measured mobilities using the computer program of O'Brien and White 4 . 

Dielectric Spectroscopy 

The dielectric and conductivity spectra were measured with the dielectric spectrometer of Myers 
and Saville 18 operated at a reference frequency of 500 Hz. The instrument calibration and operating 
procedures are described in previous publications 18 * 19 . Note that the use of a four-electrode cell 
eliminates complications due to electrode polarization 7 ®. The spectrometer's electronic drift 
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compensation system effectively controls the temperature to within 0.0001°C 18 . The absolute 
temperature inside the cell was maintained at 25.0 ± 0.1°C by a control system which employs thin film 
heaters (Minco, Minneapolis, MN) and a refrigerated air stream 21 . The low-frequency conductivity, 
which is needed in the data analysis 18 , was measured using an Altex conductivity bridge (Beckman) 
operated at 85 Hz. 

Dielectric constant and conductivity spectra for each suspension were determined at three volume 
fractions for the latex and four volume fractions in the silica experiments. After measuring the electrical 
response at a given volume fraction, a small amount of suspension was withdrawn from the cell and saved 
for a gravimetric determination. The colloidal material removed was replaced with the same volume of 
pure electrolyte yielding a suspension of lower particle concentration. Densities for polystyrene and 
silica were taken to be 1.057 and 1.8, respectively. Typical particle concentrations were 2-6% by volume. 
The reproducibility of the gravimetric determination was ± 1.0 x 10 -4 volume fraction units. 

RESULTS AND DISCUSSION 

Dielectric Spectroscopy 

The dielectric constant £ and the specific conductivity O were measured over the frequency range 0.5- 
200 kHz. The dielectric spectrometer is designed to measure the electrical properties as normalized 
frequency-difference functions, or 

[e(f,(j)) - £,]/<() , [CT(f,<t>) -CT r ]/<t> . (1) 

Here f is the linear frequency in Hz, <J> is the particle volume fraction, and the subscript r identifies 
quantities measured at the reference frequency, i.e. O r = 0(f r ,<{)). The frequency-difference functions 
represent the change in the electrical properties relative to their values at the reference frequency, 
divided by the particle volume fraction. By subtracting the signals at the reference frequency, the 
frequency-difference functions eliminate the frequency-independent contribution of the background 
electrolyte, thereby increasing the signal-to-noise ratio of the measurement 18 . For dilute suspensions. 
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the normalized frequency-difference functions are independent of <j> to the extent that particle-particle 
interactions are negligible. 

We begin by describing the behavior of a representative system in detail, to illustrate the 
relationship between the specific electrical properties and the related normalized frequency-difference 
functions. Figure la shows the dielectric constant (or relative permittivity) for the amphoteric latex in 
0.2mM HNO 3 a * three volume fractions, along with the pure electrolyte (<}> = 0) spectrum. At frequencies 
below 2kHz, £ varies slowly with frequency, as the field-induced polarization remains in phase with 
the oscillating electric field. In the frequency range 2-20 kHz, £ drops sharply, a phenomenon known as 
dielectric relaxation. This behavior is due to the relaxation time associated with the polarization 
response, which is comparable to the period of oscillation in this intermediate frequency range. At high 
frequencies, the polarization is unable to keep pace with the applied field, so the suspension dielectric 
constant approaches the value characteristic of the pure electrolyte (about 78 at 25°C). 

The normalized frequency-difference of the dielectric constant is displayed in Fig. lb. The 
normalized difference spectra superimpose, indicating that the change in the dielectric constant relative 
to the reference frequency is a linear function of particle volume fraction 18,1 ^. This implies that 
particle-particle interactions do not significantly influence the dielectric response, so the results can be 
interpreted using the existing dilute suspension theory. 

The behavior of the suspension conductivity is shown in Fig. 2a. Since all experiments were 
performed at 25°C, the separation between the data at different volume fractions is due to the electrical 
characteristics of the particles, rather than changes in the absolute temperature of the suspensions. It is 
clear from the data that the conductivity, in contrast to the dielectric constant, is a relatively weak 
function of frequency. Figure 2b shows the normalized conductivity difference spectra. Once again the 
data sets superimpose, with the exception of the results at 200 kHz. The deviations observed at high 
frequencies are systematic and appear to reflect an electronic limitation of the current spectrometer 
configuration. Consequently, in the results for latex which follow we report the conductivity spectra only 
up to 50 kHz. 



Figures 3 and 4 summarize the electrical response for the amphoteric latex in HC1 and HNO3 
solutions. The spectra are averages of measurements made at three volume fractions as described above 
and the standard deviation associated with each data point is typically 1-5%. Here, the dominant 
feature is the sensitivity to changes in the electrolyte concentration. The normalized frequency- 
difference spectra increase by nearly an order of magnitude as the ionic strength is elevated from O.lmM 
to ImM. The data for HCI and HNO3 suspensions are remarkably similar, especially at the lower ionic 
strengths. Only at ImM ionic strength are significant differences in the relaxation magnitude evident. 
Calculation of the dielectric loss from the measured conductivity differences 11,12 shows that the 
relaxation is centered between 5-7 kHz, and that the central frequency is essentially independent of 
electrolyte concentration. 

Figures 5 and 6 are the normalized frequency-difference spectra for the silica/HCl dispersions. The 
relaxation is much smaller than that observed in the latex experiments. For example, the magnitude of 
the normalized dielectric difference at 200 kHz is only 600 for silica in 0.2mM HCI, compared to 11,000 
for the latex under the same conditions. This means that the dielectric constant of a 1% silica dispersion 
changes by only 6 units over the entire frequency range. This decrease in relaxation magnitude increases 
the relative error in the frequency-difference functions, making it more difficult to measure the response 
accurately. The standard deviation of the frequency-difference data is 4-15% of the reported averages. 
Errors in the conductivity were larger at frequencies above 20 kHz. The signal-to-noise problem was 
especially severe in the dielectric measurements for 0.5mM HCl/silica at low frequencies. The small 
change in the dielectric constant below 7 kHz for this sample could not be resolved with sufficient 
accuracy; hence the data in Fig. 5 spans the frequency range 7-200 kHz. 

These results show that particles of comparable size in the same suspension medium can exhibit 
completely different electrical behavior. According to the classical electrokinetic theory, measurements 
of the ^-potential should specify the electrical response of a dispersion given the particle size and the 
ionic composition of the bathing electrolyte. In the next section, the relationship between the dielectric 
relaxation data and the measured mobilities (and ^-potentials) is discussed. 
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Electrophoresis 

Since the surface charge of amphoteric particles depends on [H + ] / it is possible that part of the ionic 
strength sensitivity revealed in the dielectric measurements results from an increase in the surface 
charge. Unfortunately, it is difficult to deduce the behavior of the surface charge from the 
electrophoretic mobility data (Table I). The mobility in O.lmM HNO3 slightly exceeds the theoretical 
mobility maximum 4 , so the potential at that maximum was used to estimate The measured mobilities 
at higher ionic strength fall short of the maxima, however the mobility-^ function is quite flat for these 
systems 4 so there is considerable uncertainty ( ca . 15%) in the value of £. Even if £ were known 
unambiguously, the relationship between £ and the surface charge needs to be specified. In general the 
electrokinetic charge differs from the surface charge due to ion localization inside the shear surface 22 . 
The structure of the region inside the shear surface is not well understood, and in fact it has been 
suggested that the location of the shear surface in latex systems may depend on the ionic strength 23 * 24 . 
In view of these difficulties, it is not possible to state whether the surface charge is constant over the pH 
range of the measurements. Nevertheless, the contrast between the relatively constant ^-potentials and 
the large increases in the frequency-difference spectra with ionic strength is quite striking. 

Silica is negatively charged for 3.3 < pH < 4; the point of zero charge is nominally at pH ~ 2.5 22 . 
Both the decrease in charge and the compression of the double layer as the pH is lowered tend to suppress 
the mobility, as demonstrated by the data in Table I. Here, in contrast to the latex experiments, the 
behavior of the ^-potential appears to be consistent with the behavior of the specific electrical 
properties: both £ and the frequency-difference spectra decrease as the pH is lowered. 

In the following section, the measured ^-potentials will be used to compare the dielectric relaxation 
data with predictions from the classical electrokinetic theory. 

COMPARISON BETWEEN THEORY AND EXPERIMENT 

Information obtained from electrokinetic experiments is generally interpreted using a model which 
specifies the relationship between a measurable quantity and an intrinsic property of the colloidal 
dispersion. For example, measurements of the electrophoretic mobility are used to estimate the 
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electrokinetic charge or ^-potential derived from the classical electrokinetic theory. The interpretation 
of bulk measurements is more complicated since the behavior results from the influence of a large number 
of particles and their associated double layers. If the suspension is not too concentrated and the ionic 
strength is not too low, then the double layers of neighboring particles will not interact significantly and 
the suspension can be treated as an ensemble of isolated particles. The model used in the electrophoretic 
mobility problem can be adapted to calculate the bulk suspension electrical properties given a suitable 
procedure for averaging local current densities and electric fields. 

Here, we summarize the model used to describe the response of an isolated spherical particle 
immersed in a dilute electrolyte in the presence of an externally applied electric field 4 ' 6 - 11,25 . For 
alternating applied currents, it is assumed that the time dependence can be represented by sinusoidal 
functions. The essential problem is to compute the field-induced perturbations in the equilibrium 
potential and ion distributions. If the external field is small in comparison to the electric field present in 
the equilibrium double layer, the problem can be attacked by constructing perturbation expansions in the 
dimensionless applied field. The basic equations are Poisson's equation for the distribution of electrical 
potential, the Stokes equations with an electrical body force for the low Reynolds number incompressible 
fluid motion, and conservation equations for the transport of mass and ionic species. The constitutive 
relations for the ion fluxes include the effects of convection, diffusion, and electromigration. The 
specification of boundary conditions completes the mathematical description of the problem. Far away 
from the charged surface, the fluid velocity tends to zero, the electric field is uniform and equal to the 
applied field, and the ion concentrations approach their bulk values. At the surface, the no-slip 
condition is assumed, the charge is taken to be constant, and the normal flux of ions into the particle is 
identically zero. 

The electrokinetic equations have been solved numerically by O'Brien and White 4 for the single 
particle electrophoretic mobility and by DeLacey and White 11 to predict the electrical response of a 
dilute suspension. These are the most comprehensive treatments of the classical theory available, and 
are not limited to thin double layers or small surface potentials. Accordingly, results from the O'Brien- 
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White and DeLacey-White methodologies will be used in the comparisons between theory and 
experiment. 

The procedure for comparing the experimental data to the theoretical predictions is as follows. The 
measured electrophoretic mobility is converted into its equivalent ([-potential using the computer 
program of O’Brien and White. Using this ([-potential in the DeLacey and White computer program 
yields predictions for the bulk suspension electrical properties as a function of frequency. The results of 
the computation are the coefficients in the density expansions 


e(f,<()) = e(f,0) + Ae’(f) <t> + ... 


a(f,<|>) = CT(f,0) + AK eff (f)<t> + ... 


( 2 ) 


correct to 0(<j>). Here, £(f,0) and G(f,0) are the specific electrical properties of the pure electrolyte. The 
functions Ae'and AK eff represent the increment in the electrical properties due to the particles and are 
directly related to the measured frequency-difference functions: 


[e(f,<|>)-e r ]/<j> = Ae’(f)-A8'(f r ) 
[cm - o r ]/(}> = AK (f) - AK (f r ) . 


(3) 


The normalized frequency-difference functions are equal to the frequency-difference of DeLacey and 
White’s A-functions. Thus, measurements of the normalized frequency-difference spectra are referenced 
to calculated frequency-difference functions based on ([-potentials obtained from electrophoresis. This 
procedure provides a critical test of the electrokinetic theory, since the theories for the mobility and the 
electrical response are self-consistent. 

Figure 7 illustrates the comparison for amphoteric latex/HNC>3 dispersions at two ionic strengths. 
Theoretical predictions using ([ = 111 mV for 0.2mM HNO3 and ([ = 98 mV for 1 mM HNO3 are displayed 



with data taken from Fig. 3. It is obvious that there is a substantial disagreement between the measured 
and calculated dielectric difference spectra. The magnitude of the relaxation response is much larger 
than predicted, especially at the highest ionic strength. Moreover, the relaxation is centered at 
frequencies well below the 30-40 kHz range predicted by the theory. The conductivity-difference data 
exhibits similar discrepancies between theory and experiment. 

It may be shown that these discrepancies cannot be eliminated by increasing the value of the 
potential or the effective particle size 13 . For example, increasing the ^-potential to a value of 250 mV 
for latex in 1 mM HNO3 yields an estimate of 9,000 for the magnitude of the normalized dielectric 
difference at 200 kHz, far less than the value of 44,000 observed experimentally. Such a large value for £ 
is exceedingly improbable. Furthermore, treating £ as an adjustable parameter does not address the 
discrepancy in the relaxation frequency, which is another distinct feature of the experimental data. 

Similar discrepancies have been observed in previous studies with latices^* 12,13 synthesized by a 
variety of surfactant-free emulsion polymerization recipes. Several carefully calibrated measurement 
techniques have been applied, and it is unlikely that all of these approaches could contain systematic 
errors large enough to explain the discrepancies. Moreover, the theory has been checked by an 
approximate analytical solution valid in the limit of small ^-potentials (< 50 mV), and the agreement 
is satisfactory 25 . In the absence of an obvious problem with either the experiments or the theory, the 
results rekindle questions about the suitability of polymer colloids as model systems for studies of 
electrokinetic phenomena. Given the dearth of published data on the dielectric relaxation of inorganic 
colloids, it seemed prudent to begin to investigate whether the discrepancies observed with latices occur 
with other colloidal systems as well. 

The comparison for the silica dispersions is shown in Fig. 8. The relatively small relaxation noted 
in Figs. 5 and 6 results in much improved agreement between the measured values and the theoretical 
predictions. For the 0.5mM HC1 /silica data set, agreement between experiment and theory is rather good 
at high frequencies. Discrepancies become apparent as the ionic strength is lowered, coincident with 
increases in the measured ^-potentials. Still, the experiment and theory are of the same order of 
magnitude, so that relatively small increases in the ^-potentials will improve the agreement. Analysis 



of the conductivity-difference data shows that there remains a discrepancy in the characteristic 
frequency, which is centered at ~5 kHz in the experiments, and between 10-20 kHz according to the 
theory. The origin of this discrepancy is not understood, and remains a subject for future study. 

CONCLUSIONS 

The results of this investigation confirm earlier studies^* which identified large discrepancies 
between theory and experiment for polymer latices. The electrical response of amphoteric latex is shown 
to be a sensitive function of the electrolyte concentration, and is virtually independent of the counterion 
for HC1 and HNO3. The strong concentration dependence of £ is not reflected in the electrophoretic 
mobilities, which depend only moderately on the ionic strength. Furthermore, dielectric spectra 
calculated using ^-potentials inferred from electrophoretic mobility disagree with the measured 
frequency-difference spectra. The theory does better with silica, which exhibits only a small 
enhancement in the electrical properties for the dilute dispersions at low pH. More experiments at 
higher pH and in different electrolytes are needed to ascertain whether this behavior is typical for 
silica. Dielectric spectroscopy appears to be sensitive to factors other than the particle size and the ionic 
strength, as demonstrated by the disparate behavior of latex and silica under the same conditions. 
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TABLE I 


Electrophoretic mobilities and ^-potentials for amphoteric latex and silica 
dispersions as a function of electrolyte concentration. 


Amphoteric Latex: 


[HN 03 ] 

aK 

mobility 3 

^-potential* 5 

.0001 M 

3.7 

4.5 

144 c 

.0002 M 

5.3 

4.3 

111 

.0005 M 

8.3 

4.6 

125 

.001 M 

11.8 

4.7 

98 

Silica: 

[HC1] 

aK 

mobility 

^-potential 

.0001 M 

5.4 

-2.4 

-42 

.0002 M 

7.7 

-2.3 

-37 

.0005 M 

12.1 

-1.7 

-27 


a Mobility units are pm-cm/ volt-sec. 

^-potential units are mV. 

c The ^-potential at the mobility maximum, which is 4.31 micron-cm/volt sec. 
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FIGURE CAPTIONS 


Figure la: Dielectric constants for an amphoteric latex in 0.2mM HNO3 at 25°C. 

• <{> = 0.000; ♦ <|> = 0.016; A 4> = 0.023; ■ 4> = 0.032. 

Figure lb: Normalized dielectric difference spectra for an amphoteric latex in 0.2mM HNO3 at 25°C. 

• $ = 0.016; A $ = 0.023; ■ <J> = 0.032. 

Figure 2a: Specific conductivities for an amphoteric latex in 0.2mM HNO3 at 25°C. 

• 4. = 0.000; ♦ <j> = 0.016; A <}> = 0.023; ■ <> = 0.032. 

Figure 2b: Normalized conductivity difference spectra for an amphoteric latex in 0.2mM HNO3 
at 25°C. ♦ $ = 0.016; A <j> = 0.023; ■ <}> = 0.032. 

Figure 3: Normalized dielectric difference spectra for an amphoteric latex as a function of ionic 

strength at 25°C. Data points are averages of measurements made at three volume fractions. 

• O.lmM HNO3; ■ 0.2mM HNO3; A 0.5mM HNO3; ♦ ImM HNO3; O O.lmM HC1; 

□ 0.2mM HC1; A 0.5mM HC1;0 ImM HC1. 

Figure 4: Normalized conductivity difference spectra for an amphoteric latex as a function of ionic 

strength at 25°C. Data points are averages of measurements made at three volume fractions. 

• O.lmM HNO3; ■ 0.2mM HNO3; A 0.5mM HNO3; ♦ ImM HNO3; O O.lmM HC1; 

□ 0.2mM HQ; A 0.5mM HC1; O ImM HC1. 

Figure 5: Normalized dielectric difference spectra for silica as a function of ionic strength at 25°C. 

Data points are averages of measurements made at four volume fractions. 

■ O.lmM HC1; A 0.2mM HC1; ♦ 0.5mM HC1. 

Figure 6: Normalized conductivity difference spectra for silica as a function of ionic strength at 25°C. 

Data points are averages of measurements made at four volume fractions. 

■ O.lmM HC1; A 0.2mM HC1; ♦ 0.5mM HC1. 

Figure 7: Comparison between experiment and theory for an amphoteric latex in HNO3. ' 

ImM HNO3: A experiment, theory; 0.2mM HNO3: ■ experiment, theory. 

Figure 8: Comparison between experiment and theory for silica in HC1. 

O.lmM HC1: ■ experiment, theory; 0.2mM HC1: A experiment, theory; 

0.5mM HC1: ♦ experiment, theory. 


1 5 




500 - 



CO 





F-COlOtFCOCVI't-iO 
00 CO GO CD 00 GO CO CD 

(rno/sr/) d 






Figure 


o ♦ 04 no 

O* < BO 

O ♦ <1 

\ 

<o+> mm* 




\~mm 

O 

['TTTT |'T I 1 1 1 

o 

II 1 1 1 | 1 1 II ! 

O 

i i i ii| i in i 
O 

i ii i i | i mi 
O 

i i i i if 

o 

LO 

o 

in 

o 

in 


CV2 

CV2 

rH 





(rno/gr/) d)/[*n - (^‘j)d] 




[c(f ,<p) - <T r ]/(p (fjS/ C 


Figure 6 



l 


Frequency(kHz) 



- Zt)\/cP 


Figure 7 






The Dielectric Response of Polystyrene Latexes: 
Effects of Alterations in the Structure of the Particle Surface 

Manuscript #8243K-revised March 28, 1990 

■ o cx p petoy )ys 2 • Coltoic* 


L. A. Rosen and D. A. Saville 

Department of Chemical Engineering 
Princeton University 


Princeton NJ 08544 




Abs tra ct 

or t'©t>3V\»C£5' 

The presence of a "hairy" layer- has been hypothesized to explain 
discrepancies between the electrokinetic theory and experiments with polystyrene 
latexes by Van der Put and Bijsterbosch [J. Colloid Interface Sci. 92, 499 (1983)]. 
We tested the hairy4ayei^hypothesis by making electrokinetic measurements on 
latexes before and after they had been heated above the polystyrene glass transition 

Svncxo’H'i +v>e pasr+iile, b ^ aUoivi'ocj 

temperature. This sort of treatment should^induce the charged ■ polymer filaments 
which comprise the layer to collapse onto the particle surface. Low-frequency 
dielectric measurements reveal that heat-treatment improves the agreement 
between experiment and theory. For an amphoteric latex, the dielectric constant 
and conductivity decrease with heat-treatment time, and the characteristic 
relaxation frequency increases. This is consistent with the hairy surface 
hypothesis, suggesting that heat-treatment smooths the surface, creating particles 
which more closely conform to the assumptions of the classical theory. 
Experiments with an anionic latex further suggest the presence of surface 
structure, but discrepancies between experiment and theory persist; it appears 
that factors other than hairiness contribute to these discrepancies. The 
suspension’s dielectric properties are evidently very sensitive to the state of the 
particle surface, and thus dielectric spectroscopy offers a promising means of 
studying the structure of colloidal interfaces. 
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Introduction 


Polystyrene latexes produced by surfactant- free emulsion polymerization (1) 
consist of monodisperse spherical particles stabilized by covalently-bound charge. 
These dispersions should be excellent model systems for studies of colloid stability 
and electrokinetic phenomena. Nevertheless, a large body of experimental 
evidence indicates that latexes do not behave as the standard or ’’classical" 
theories predict. For example, results of coagulation studies disagree with the 
theoretical stability ratio-ionic strength relationships (2,3). The ionic strength 
dependence of the ^-potential for certain latex systems exhibits a maximum not 
predicted by electrokinetic theory (4-6). Conductivity measurements on polystyrene 
latex porous plugs provide ^-potentials which are inconsistent with ^-potentials 
determined from microelectrophoresis (7,8). Furthermore, the ^-potentials 
obtained from microelectrophoresis differ from those inferred from dilute 
suspension conductivity measurements (8,9). 

Pervasive inconsistencies are also apparent in the dielectric properties of 
polystyrene latexes at low frequencies. Here, the measured dielectric constant 
can exceed the theoretical prediction by an order of magnitude or more, and the 
relaxation frequency can be much lower than predicted (10). 

The source of the disagreement between experiment and theory is not well 
understood, although it appears likely that the discrepancies arise from effects 
associated with the particle surface rather than the diffuse double layer. The 
classical theories are based on the Gouy-Chapman theory of the equilibrium 
diffuse layer, as described by the full non-linear Poisson-Boltzmann equation. The 
available experimental evidence indicates that the Gouy-Chapman theory 
quantitatively predicts the electrostatic repulsion between two charged mica 
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surfaces, down to separations of a few nanometers (11,12). Thus it is believed that 
this theory provides an adequate representation of the diffuse layer. 

Studies of electrode and metal oxide surfaces show that there is a region close 
to the surface where the ionic composition and potential distribution differ from 
that of the diffuse layer. This region, called the Stem layer, consists of specifically 
adsorbed ions and solvent molecules (13). In the classical theories, the charge 
inside the Stern layer is constant and the Stern layer ions are immobile. 
Accordingly, the layer is static, and fully decoupled from transport processes 
occuring in the diffuse layer. It was shown recently that accounting for non- 
equilibrium effects in the Stern layer can significantly affect the predicted 
electrokinetic behavior of colloidal suspensions (14,15). 

A related aspect called into question by recent experiments (6,16) is the 
smoothness of the polystyrene surface. It has been suggested that the surface may 
be "hairy" due to the presence of terminally-charged polymer chains which 
interact electrostatically at low to moderate ionic strengths (6). Should such 
chains protrude into the Stern layer or the inner diffuse layer, electrokinetic effects 
will be altered in several ways. The displacement of covalently-bound charge from 
the surface means that the distribution of fixed charge will differ from that 
specified in the classical theory which assumes that the fixed charge is located in 
a region of infinitesimal thickness at the particle surface. Extended chains will 
shift the location of the shear surface away from the interface, resulting in a 
decrease in the ^-potential relative to the value for a smooth surface possessing the 
same fixed charge. This expansion of the inner double layer should accentuate 
the importance of the dynamic Stern layer effects (14,15) not accounted for in the 
classical electrokinetic theory. Interfacial ion transport will be hindered by 
surface hairs, so that surface ion mobilities will be lower than their corresponding 
bulk values. 
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Suppose we adopt the hypothesis that the polystyrene surface is hairy. If this 
is true, then the surface should be susceptible to heat-treatment above the 
polystyrene glass transition (16). Under these conditions, the extended chains will 
be mobile, and may be expected to collapse onto the particle surface. Such 
alterations in the structure of the interfacial region should be revealed in the 
electrokinetic behavior, for the reasons discussed above. 

We have made low-frequency dielectric spectroscopy and 
microelectrophoresis measurements on "untreated" and "heat-treated" 
polystyrene latexes. The untreated latexes exhibit large deviations from the 
predicted response, but heat-treatment improves agreement with the predictions of 
the classical theory. With an amphoteric latex, four hours of treatment produces a 
suspension which behaves approximately as the classical theory predicts. These 
results are consistent with the hairy surface hypothesis. With an anionic latex, 
heat-treatment for twelve hours improves the agreement between experiment and 
theory, but significant discrepancies remain. 

The presentation is organized as follows. In the next section, we review the 
basic features of the classical or standard theory of electrokinetics with which our 
results will be interpreted. Following this is a description of the experimental 
protocols. Next we outline the techniques used to depict the results, discuss the 
use of difference spectra, and explain how experiment and theory will be 
compared. The results are described next, followed by a set of concluding 
remarks. 
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Theoretical Considerations 


In what follows, the terminology classical theory or standard electrokinetic 
model means the set of partial differential equations and associated boundary 
conditions used to describe the behavior of an isolated spherical particle i mm ersed 
in a dilute electrolyte (17,18), in the presence of an externally applied electric field. 
The equations are Poisson's equation for the distribution of electrical potential, the 
Stokes equations with an electrical body force for the low Reynolds number 
incompressible fluid motion, and conservation equations for the transport of mass 
and ionic species. The constitutive relation for the ion transport includes the 
effects of convection, diffusion, and electromigration. Far away from the surface, 
the influence of the particle disappears, and the field quantities obtain uniform 
values. At the surface, the no-slip condition is assumed, the charge is taken to be 
constant, and the normal flux of ions into the particle is identically zero. Note that 
this formulation does not take into account effects such as surface charge 
polarization and Stem layer transport (14,15). 

The standard electrokinetic model equations have been solved numerically by 
O’Brien and White (19) for the single particle electrophoretic mobility and by 
DeLacey and White (20) to predict the dielectric response of a dilute suspension. 
These are the most comprehensive treatments available, and are not restricted to 
thin double layers or small surface potentials. For our suspensions, which 
consist of small, highly charged particles at low ionic strengths, more 
approximate models (21-24) are not generally applicable. Accordingly, results 
from the O’Brien-White and DeLacey-White methodologies will be used in the 
comparisons between theory and experiment. 
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Experimenta l Met h ods 


An amphoteric polystyrene latex, designated as Latex A, was synthesized 
according to the procedure of Homola and James (25), using the modification 
discussed by Harding and Healy (26). Based on the synthesis recipe, the predicted 
isoelectric pH for this latex is approximately 5.4. For our experiments, the latex 
was in the cationic form. 

The anionic latex, Latex B, was synthesized via the method of Goodwin et al. 
(27). The electrophoretic mobilities of Latex B were uniformly negative, consistent 
with the hypothesis that the charge arises from the dissociation of surface sulfate 
groups and possibly some carboxylic acid groups. 

The latexes were cleaned by sedimentation/decantation in a Beckman Model 
L5-65 ultracentrifuge with a Type 21 rotor. The first four cycles consisted of an acid 
pre-wash (ImM HC1). Subsequent wash cycles employed the electrolyte to be used 
as the suspension medium; normally ten cycles were required for the conductivity 
of the supernatant to reach a value within 1% of the wash electrolyte. The 
untreated Latex A was easily resuspended without sonication. Heat-treated Latex 
A was more difficult to resuspend, requiring light sonication in a bath-type 
ultrasonic cleaner (L & R Model T-9B, Sargent Welch). Latex B required sonication 
both before and after heat-treatment. 

Before heat-treatment the latex was diluted to a volume fraction of 0.5% solids, 

i 

and placed in a citrate glass beverage bottle with a teflon-lined bottle cap. The 
purpose of the dilution was to minimize aggregation during the high-temperature 
treatment. Such aggregates would effectively increase the polydispersity, making 
the dielectric data difficult to interpret. The bottled suspensions were heated to 120 
± 2 °C in a Deluxe pressure cooker/canner (Presto Ind., Eau Claire WI) equipped 
with a dial thermometer (Omega, Stamford CT). At the end of the treatment, the 
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suspension was allowed to cool slowly to room temperature inside the pressure 
cooker. The heat-treated suspensions showed no evidence of particle coagulum. 
We did observe that rapid cooling, e.g. immediate refrigeration, caused a number 
of large fused aggregates to form. 

After heat-treatment, the latex was reconcentrated (2-6% solids by volume) by 
centrifugation, with additional cycles as needed to bring the conductivity to within 
1% of the value for the pure electrolyte. The last cycle always took place between 12 
and 24 hours before the dielectric spectroscopy experiments. 

Particle size and monodispersity were determined by transmission electron 
microscopy. The particles appeared uniform and spherical before and after heat- 
treatment. The magnification was calibrated with latex standards of known 
diameter (Interfacial Dynamics, Portland OR). The results (Table I) indicate that 
the mean particle diameter and uniformity were essentially unaffected by the heat- 
treatment procedure. For Latex B, a few large spherical particles were observed in 
the electron micrographs following heat-treatment. These appear to have been 
formed by the coalescence of primary particles above T g . 

Electrophoretic mobilities were measured with d Rank Brothers Mark II 
instrument, in the four-electrode cell configuration with palladium current 
electrodes. Mobilities were converted into ^-potentials using the computer program 
of O’Brien and White (19). 

The dielectric and conductivity spectra were determined using the four- 
electrode frequency-difference dielectric spectrometer of Myers and Saville (28) 
operated at a reference frequency of 500 Hz. The instrument operation and 
calibration procedures are described in detail elsewhere (28,29). Note that the four- 
electrode cell eliminates complications due to electrode polarization (30) and the 
electronic drift compensation system provides precise "effective" temperature 
control, estimated to be better than 0.0001°C (28). The absolute temperature inside 
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the cell was maintained at 25.0 ± 0.1°C by a temperature control system which 
employs thin film heaters (Minco, Minneapolis, MN) and a refrigerated air stream 
(31). The low-frequency conductivity, which is needed in the data analysis (28), was 
measured using an Altex conductivity bridge (Beckman) operated at 85 Hz. 

The normalized difference spectra for each suspension were determined at 
three or four volume fractions using a successive dilution technique. After 
measuring the dielectric response at a given volume fraction, a small amount of 
suspension was withdrawn from the cell and saved for a dry-weight determination, 
assuming the density of polystyrene to be 1.057. The material removed was 
replaced with the same volume of pure electrolyte yielding a suspension of lower 
particle concentration. 
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Methods of Presenting the Experimental Data 


Dielectric and conductivity data are presented in terms of normalized 
frequency-difference functions (28), 


[e(f,<|>) -£,]/<!> , [o(f,<|>) -c r ]/<> 


[ 1 ] 


where e r and a r are the dielectric constant and conductivity at the reference 
frequency, i.e. e(f r ,<J>) and o(f r ,4>), and <j> is the particle volume fraction. By 
subtracting the signals at the reference frequency, the frequency-difference 
functions eliminate the frequency-independent contribution of the background 
electrolyte. The normalized frequency-difference functions are independent of <}) to 
the extent that particle-particle interactions are negligible. These functions are 
directly related to DeLacey and White's (20) functions, which are the O (<}>) 
coefficients in the density expansions: 


e(f,4>) = e(f,0) + Ae’(f) § + ... 

CF(f,<(>) = o(f,0) + AK eff (0<ti + ... . 


[ 2 ] 


Here, e(f,0) and o(f,0) refer to the dielectric constant and conductivity of the pure 
electrolyte (<{> = 0), which are frequency-independent over the measurement range 
(18). It follows that 


[e(f,4>) - £r]/<{> = Ae'(f) - Ae'(f r ) 
[o(f,0) - Or ]/<| ) = AK eff (0 - AK eff (f r ) . 


[3] 
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The normalized frequency-difference functions employed in this work are 
equivalent to the frequency-difference of the DeLacey-White functions. 

Dielectric relaxation data is sometimes presented in terms of a complex 
dielectric constant e’ + ie H . Here, e’ is the real or "in-phase" part of the dielectric 
constant; e" is the imaginary or "out-of-phase" component, referred to as the 
dielectric loss. In a perfect dielectric e" is zero; nonzero values arise in colloidal 
systems due to diffusion, which reduces the orientation polarization and therefore 
the dielectric constant e' (18). The dielectric loss can be related to the measurable 
suspension conductivity: 


o(f,4>) =o(0,<{))+ co e v e "(f,<j)). [4] 

Here, co represents the angular frequency in radians/second; f = co/2tt. Equation [4] 
separates the conductivity into a dc component o(0,<{0 and a frequency-dependent 
term which involves the dielectric loss function. The constant e v denotes the 
vacuum permittivity, 8.854 x lO-^F/m. At the reference frequency, Eq. [4] can be 
rewritten: 

o r (f r ,4>) =o(0,<}))+ C 0 r e v e r "(fr>40 . [5] 

Subtracting [5] from [4] eliminates the dc conductivity, yielding an expression for 
the conductivity difference function: 


o(f,4>) - o r (f r ,4>) = co e v e M (f,4>) - co* e v e r "(f r ,4>) . 


[ 6 ] 



Now if the conductivity measured at the reference frequency o r (f r ,<j>) is equal to the 
dc conductivity o(0,<|)), then it follows from Eq. [5] that ^"(fr,^) = 0. According to Eq. 
[ 6 ], 


e"(f, <»/<() = [c(f,<{>) - 0 r (f r ,<t>)]/co e v . [7] 

Measurements show that the suspension conductivity is essentially constant below 
1 kHz, so the assumption a r (f r ,<{)) = c(0,<{>) is justified. The dielectric loss function 
exhibits peaks in the vicinity of the characteristic relaxation frequencies, and thus 
provides a convenient means of viewing changes in the relaxation response. 

The frequency-difference functions for each experiment were checked to 
establish that the data were linear with particle volume fraction. Correlation 
coefficients from linear regression normally exceeded 0.99. The conductivity- 
differences above 50 kHz were non-linear in <|> so we report conductivity data only up 
to 50 kHz. For each amphoteric suspension, the reported normalized difference 
functions reflect average values of the successive dilution experiments at three 
volume fractions. Measurements were made at four volume fractions in 
experiments with anionic latex. Typically, the standard deviation was 1-5% of the 
mean. The only exception was with the four hour heat-treated amphoteric latex. 
In this case, the relative error was larger due to the small magnitude of the 
signals, as will be discussed below. 



Experimental Results 
Amphoteric particles (cationic form) 

Figure 1 shows the normalized dielectric difference function for Latex A in 
ImM HC1 as a function of treatment time at 120 °C. The initial one hour treatment 
produced a remarkable change in the dielectric response in that the magnitude of 
the high frequency dielectric difference decreased from almost 49,000 for the 
untreated latex to 17,000. The normalized conductivity function, shown in Fig. 2, 
underwent a similar depression. 

This data must be viewed in conjunction with the electrophoretic mobility, 
which also exhibited interesting behavior. The mobility (Fig. 3) increased from 5.13 
to 6.18 jim-cm/volt-sec after one hour, which substantially exceeds the theoretical 
maximum predicted by the O'Brien- White solution (19). It is difficult to interpret 
this data in terms of ^-potentials, since, according to the theory, there is no £- 
potential which corresponds to the measured mobility. Critical coagulation 
concentrations (CCC) were also measured for the latexes. For the untreated Latex 
A, 0.25M < CCC < 0.30M in a KNO3/HNO3 solution at pH=1.3. The one hour heat- 
treated Latex A had lower aggregation stability, with .20M < CCC < .25M at pH=1.3. 
So, despite the increase in mobility, the heat-treated latex was more susceptible to 
aggregation. 

Additional heat-treatments resulted in further attenuation of the low 
frequency relaxation signals. Figure 4 shows the behavior of e7<|> as a function of 
treatment time. Values of the loss function were computed from the measured 
conductivity differences using Eq. [7]. From Fig. 4, it is evident that the low 
frequency response (<10 kHz) is supressed by the heat-treatment procedure. At the 
same time, there is a suggestion that a higher frequency mode (-50 kHz) is present. 
After four hours of treatment, there is no longer a low-frequency peak, and the 
principal relaxation occurs above 20 kHz. The small magnitude of the conductivity- 
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difference signal combined with the inherent limitations of the measurement 
technique at high frequencies (29,30) made it impossible to resolve the high 
frequency peak for the four hour data. 

Returning to Fig. 1, we see that the magnitude of dielectric-difference after 
four hours at 120°C is only 2300. The standard deviation associated with the value 
of [e(f,4>) - e r ]/<{> at 200 kHz is about 15% of the mean. This relatively large error is 
due to the decrease in the dielectric constant of the suspension brought about by 
four hours of heat- treatment. For a suspension containing 2% particles by volume, 
the dielectric constant at 500 Hz is only about 125, compared to over 1000 for the 
untreated latex. This huge decrease in signal means that it is much more difficult 
to measure the dielectric-difference function (29,30). Nevertheless, we are able to 
achieve sufficient accuracy to make meaningful comparisons with the 
electrokinetic theory. 

The electrophoretic mobility of the four hour heat-treated particles (Fig. 3) 
decreased to within the range permitted by the O’Brien-White solution. The 
measured mobility corresponds to a ^-potential of 81 mV, which was used as input 
in the DeLacey-White computer program. The predicted response is shown in 
Figures 1-2. It is clear that the agreement between experiment and theory has 
been vastly improved by the heat-treatment procedure. The measured and 
predicted high frequency magnitudes agree to within the experimental error (15%), 
and the relaxation frequency is much closer to the value expected from the classical 
theory. Thus, the electrokinetic theory appears internally consistent when applied 
to amphoteric particles heat-treated for a sufficiently long period of time. 

Based on the mobility measurements, there is considerable uncertainty in the 
^-potential as a function of treatment time. In particular, the ^-potential for 
untreated and one hour heat-treated latexes may be higher than the 81 mV 
estimate for the four hour latex. A sensitivity analysis (28,29) demonstrates that 



even after increasing the ^-potential to 250 mV, the classical theory fails to predict 
the dielectric data for untreated and one hour heat-treated latexes. Thus, the large 
discrepancies noted in Figures 1,2, and 4 are not due to ambiguities in the £- 
potential. 

Anionic Particles 

Following the experiments with Latex A, an anionic latex was studied to see if 
the sensitivity to heat-treatment was a peculiarity of the amphoteric surface or a 
more general characteristic of polystyrene latexes. Figures 5-7 summarize the 
dielectric spectroscopy experiments with the anionic latex; the electrophoretic 
mobility history is displayed in Fig. 8. It is apparent that Latex B does not respond 
as rapidly to heat-treatment as Latex A. Four hours of treatment, which 
dramatically transformed Latex A, results in a modest 40% reduction in the 
magnitude of the dielectric relaxation for Latex B. Additional heating produces 
smaller shifts, with very little difference between the eight hour and twelve hour 
data. Furthermore, the shift to higher relaxation frequencies observed with Latex 
A apparently does not take place with Latex B. From Fig. 7, the low-frequency 
peak in e7<J) persists even after twelve hours of treatment, and if any shift occurs it 
is to slightly lower frequencies. 

The electrophoretic mobility (Fig. 8) increases by about 15% after heat- 

treatment to values which approach or exceed the theoretical maximum. The £- 

> 

potential at the maximum is -130 mV. Using this ^-potential, the theoretical 
spectra were computed for Latex B (Figs. 5-7). Even after twelve hours of heat- 
treatment, a significant discrepancy remains between the measured response and 
the theoretical prediction. It is clear from the behavior of Latex B that heat- 
treatment is not always sufficient to produce classical or "model" polystyrene 
particles. 
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Discussion 


The concept of a hairy layer on the polystyrene surface was advanced by Van 
der Put and Bijsterbosch (6) to explain discrepancies between theory and 
experiment observed in their electrokinetic data. According to these authors, the 
surface is covered by a flexible hairy layer consisting of polymer chains terminated 
by charged groups which derive from the initiator used in the latex synthesis. The 
thickness of the layer is determined by interactions between the charged end 
groups, and thus depends on the ionic strength. Van der Put and Bijsterbosch 
estimated a layer thickness of 1-2 nanometers at low ( £ 1 mM) ionic strength (6). 
Convection is presumed absent in the hairy layer, however ionic conduction within 
the layer makes an important contribution to the electrokinetic response (6). 

Independent non-electrokinetic evidence for the hairy layer model was 
reported by Goossens and Zembrod (32). The hydrodynamic radius of latex 
particles was measured by photon correlation spectroscopy, and was shown to be a 
sensitive function of pH and ionic strength under conditions where particle-particle 
interactions were unlikely. The hydrodynamic radius varied by as much as 12 
nanometers; at the same time electron microscopy detected no dependence of 
particle size on pH. The variation in hydrodynamic radius was attributed to 
changes in the surface layer thickness due to the pH-ionic strength dependence of 
electrostatic segment-segment interactions (32). 

Recently, Chow and Takamura (16) tested the hairy layer hypothesis by heat- 
treating a carboxylated styrene-acrylamide copolymer latex at 115 °C for 16 hours. 
Electrophoretic mobilities were measured with a laser doppler apparatus. The 
results showed that the ^-potential increased significantly following the heat- 
treatment. This observation was attributed to a shift in the shear plane due to the 
creation of a smoother surface, resulting in an increase in the electrokinetic 


15 



charge. It was concluded that the surface roughness or hairiness was an 
important factor in determining ^-potentials from electrophoretic mobility 
measurements (16). 

The present work shows how another electrokinetic property, the dielectric 
response, is influenced by the structure of the particle surface. Previous studies 
with untreated anionic (33,34) and amphoteric (10) latexes identified very large 
discrepancies between experiment and theory. Based on these results, it appeared 
that either the electrokinetic theory is flawed, or else the "model" polystyrene 
suspensions do not conform to the assumptions of the theory. Considering the 
second possibility, attention is naturally directed to the particle surface. From a 
structural standpoint, less is known and more is assumed about the surface than 
either the diffuse layer or the adjoining bulk phases. Our results suggest that the 
structure of the surface region can be altered, but the particle size remains 
constant as a result of heat-treatment. 

The results of our experiments show that the dielectric response is very 
sensitive to the state of the particle surface. Indeed, the changes in the dielectric 
properties with heat-treatment are much more pronounced than the corresponding 
shifts in the electrophoretic mobility reported in this work or in the previous heat- 
treatment study (16). The magnitude of the frequency-difference spectra for Latex 
A decreases continuously with treatment time, reaching values within 15% of the 
theoretical prediction. The decrease in the magnitude is accompanied by a shift in 

j 

the relaxation to higher frequencies. The frequency shift indicates that ion 
transport near the interface is facilitated by the change in surface structure, since 
scaling arguments predict that the characteristic frequency is directly proportional 
to the ion diffusivity (18). The apparent increase in the ion diffusivity can be 
explained in light of the hairy layer model: the collapse of surface hairs decreases 
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the resistance to surface transport, and this leads to an increase in the ion 
diffusivity, which is reflected experimentally by the shift to higher frequencies. 

The behavior of the Latex A mobility as a function of treatment time (Fig. 3) 
draws attention to the complex nature of the amphoteric latex surface. The initial 
increase in mobility to values above the theoretical maximum is similar to the 
behavior observed by Chow and Takamura (16) for the carboxylated latex. The 
subsequent decrease for longer times suggests that effects other than surface- 
smoothing may be operative. One possibility is that the heat-treatment removes 
some of the fixed charge on the particle by hydrolysis. This possibility will be 
checked in future studies by determining the titratable charge before and after 
heat-treatment. The increase in ion mobility apparently brought about by surface 
smoothing might give rise to increased electroosmotic drag on the particle, which 
would tend to lower the electrophoretic mobility. It appears likely that a number of 
surface effects contribute interactively to the electrokinetic response of amphoteric 
latex. 

The experiments with Latex B demonstrate that surface hairiness by itself 
cannot explain all of the discrepancies between theory and experiment. As with 
Latex A, the magnitude of the relaxation response decreases with treatment time, 
however much longer treatments are required to affect the signals. Inspection of 
Fig. 7 shows that the characteristic frequency, which is lower than the theoretical 
prediction by about a factor of two, does not shift to higher frequencies even after 
twelve hours of heat-treatment. The slight shift to lower frequencies may have been 
caused by the presence of a relatively small number of large particles formed 
during the heat-treatment; such particles were occasionally observed in the 
electron micrographs of heat-treated Latex B. 

The electrophoretic mobility of Latex B (Fig. 8) is relatively insensitive to the 
heat-treatment. This is reminiscent of earlier experiments by Ma et al. (35), who 


demonstrated that for latex particles of low (< 5 pC/cm 2 ) titratable charge density, 
complete hydrolysis of titratable sulfate groups to the hydroxyl form did not sensibly 
affect the electrophoretic mobility. The mobilities and ^-potentials measured for 
Latex B are considerably higher than those of Ma et al., so if hydrolysis did take 
place the surface was altered in a manner which enhanced the development of 
negative charge, presumably by chloride ion adsorption. 

Finally, it is noteworthy that the dielectric response decreases concommitant 
with the small increases in mobility for Latex B. Theory predicts that the 
magnitude of the dielectric response should increase monotonically with the 
potential, so clearly there is some difficulty in using ^-potentials from 
microelectrophoresis to predict dielectric relaxation of polystyrene latexes. This 
finding is consistent with the earlier data of Van der Put and Bijsterbosch (7) and 
Zukoski and Saville (9), demonstrating once again that the ^-potential is 
insufficient to characterize the electrokinetic behavior of polystyrene latexes. The 
use of surface structural models to account for ion binding and Stern layer 
transport (14,15) may help explain the discrepancies between theory and 
experiment in the dielectric spectra. 

Conclusions 

Dielectric measurements on polystyrene latex suspensions indicate that 
treatment above T g alters the structure of the particle surface. The effects of the 

heat-treatment increase systematically with treatment time. The results for an 
amphoteric latex are consistent with the view that the surface of these particles is 
hairy, and that the heat-treatment is able to collapse the surface hairs, creating 
particles which resemble the smooth spheres assumed in the classical theory. For 
our anionic latex, the agreement between theory and experiment is greatly 
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improved, however persistant discrepancies in the relaxation frequency suggest 
that other physicochemical factors are present. These require elucidation. Our 
results concur with earlier studies (7-10) which demonstrate that the ^-potential 
measured from microelectrophoresis does not adequately characterize the 
electrokinetic behavior of polystyrene latexes. Dielectric spectroscopy, which 
provides more information per experiment than traditional electrokinetic 
methods, also appears to be more sensitive to changes in the structure of the 
particle surface than microelectrophoresis. 

There are, no doubt, many reasons for the non-classical behavior of latex 
particles. Among these may be ion adsorption of various sorts, surface roughness 
or porosity, and the presence of polymer filaments that form a hairy surface layer. 
For example, van den Hul and Vanderhoff (36) found that heat- treatment reduced 
the porosity of latex particles based on differences between surface areas derived 
from electron microscopy and the BET method. Certainly the mix of factors will 
vary according to the synthesis procedure and subsequent treatment of the latex. 
Our results show that heat-treatment alters electrokinetic properties substantially 
and in a manner consistent with the presence of a hairy layer on the amphoteric 
particles studied. Nevertheless, much more work will be necessary to establish 
the roles of the various mechanisms with a given particle and our results do not 
rule out other explanations. Indeed, in the case of the anionic particles, non- 
classical behavior persists in the face of heat-treatment. 
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TABLE I 


Particle size information for untreated and heat-treated polystyrene latexes 



diameter(nm) 

% variation 

Latex A, untreated 

226 

3.0% 

Latex A, heat-treated 4 hr. 

229 

1.5% 

Latex B, untreated 

416 

1.3% 

Latex B, heat-treated 4 hr. 

417 

2.6% 

Latex B, heat-treated 12 hr. 

421 

1.6% 


22 




Figure 1: The effect of heat-treatment on the normalized dielectric-difference 

spectra of Latex A in ImM HC1. • untreated; ■ 1 hour at 120 °C; 

A 2 hours at 120 °C; ♦ 4 hours at 120 °C; prediction from the 

classical theory, £ = 81 mV. 


Figure 2: The effect of heat-treatment on the normalized conductivity-difference 

spectra of Latex A in ImM HC1. • untreated; ■ 1 hour at 120 °C; 

A 2 hours at 120 °C; ♦ 4 hours at 120 °C; prediction from the 

classical theory, £ = 81 mV. 


Figure 3: The effect of heat-treatment on the electrophoretic mobility of Latex A 

in 1 mM HC1. Error bars reflect one standard deviation about the 
mean. Theoretical maximum according to the O'Brien and White 
numerical solution (19). 


Figure 4: The effect of heat-treatment on the normalized dielectric loss spectra 

of Latex A in ImM HC1. • untreated; ■ 1 hour at 120 °C ; 

A 2 hours at 120 °C; ♦ 4 hours at 120 °C; prediction from the 

classical theory, £ = 81 mV. 


Figure 5: The effect of heat- treatment on the normalized dielectric-difference 

spectra of Latex B in ImM KC1. • untreated; ■ 4 hours at 120 °C; 

A 8 hours at 120 °C; ♦ 12 hours at 120 °C; prediction from 

the classical theory, £ = -130 mV. 
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Figure 6: 


Figure 7: 


Figure 8: 


The effect of heat- treatment on the normalized conductivity-difference 
spectra of Latex B in ImM KC1. • untreated; ■ 4 hours at 120 °C; 

A 8 hours at 120 °C; ♦ 12 hours at 120 °C; prediction from 

the classical theory, £ = -130 mV. 

The effect of heat- treatment on the normalized dielectric loss 
spectra of Latex B in ImM KC1. • untreated; ■ 4 hours at 120 °C; 

A 8 hours at 120 °C; ♦ 12 hours at 120 °C; prediction from 

the classical theory, £ = -130 mV. 

The effect of heat- treatment on the electrophoretic mobility of Latex B 
in 1 mM KC1. Error bars reflect one standard deviation about the 
mean. Theoretical maximum according to the O'Brien and White 
numerical solution (19). 
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Electrokinetic Properties of Particles in Concentrated Suspensions: 

Heterogeneous Systems 
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The first paper in the series [J. Colloid Interface Science 115, 422 ( 1987)] dealt with the electrokinetic 
behavior of homogeneous suspensions. There it was found that the net particle-particle interaction was 
negligible for particles undergoing electrophoresis; the particles behaved as electromechanical ghosts. In 
the current work, attention is focused on heterogeneous systems where a portion of the particles that 
make up the suspension have electrokinetic properties different from the remainder. Transparent particles 
(fast ghosts) were prepared by lysis of red blood cells to remove hemoglobin and fixed using glutaraldehyde. 
To prepare particles with a different mobility (slow ghosts), the fixed ghosts were treated with neur- 
aminidase to remove part of the sialic acid. Tracer particles were prepared in similar ways by omitting 
the lysis step. The electrical conductivities of suspensions of fixed ghosts and neuraminidase treated fixed 
ghosts were almost indistinguishable, whereas the electrophoretic mobility of the latter particles was less 
than half that of the untreated particles. Although the average conductivity of a 50-50 mixture of fast 
and slow ghosts decreased rapidly with increasing particle concentration, the mobilities of individual 
constituents did not. To within the experimental accuracy of the measurement, the mobilities of individual 
constituents varies as ( I - C), C being the volume fraction of particles, irrespective of the electrokinetic 
makeup of the suspension. This behavior is the same as that encountered with homogeneous systems. 
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I. INTRODUCTION 

Electrokinetic processes are used to char- 
acterize all sorts of suspended particles and 
play important roles in the fractionation of 
mixtures into relatively pure subpopulations 
(1,2). However, despite their importance, few 
experimental studies exist on heterogeneous 
mixtures. Most studies deal with homogeneous 
suspensions where average mobilities are 
measured by moving boundary or mass trans- 
port techniques (3-8). These yield little in- 
formation about the dynamics of individual 
constituents. 

Modeling of suspension electrokinetics has, 
in fact, preceded detailed experimental studies. 
The cell models of Moller, Van Os, and Ov- 
erbeek ( 7, 8 ) and Levine and Neale ( 9 ) deal 

1 Now at the Department of Chemical Engineering, 
University of Illinois, Urbana, IL 61801. 

2 To whom correspondence should be addressed. 


with particle concentration effects in homo- 
geneous suspensions, whereas Reed and Mor- 
rison (10) analyzed the motion of spheres with 
different zeta potentials. Anderson (11) used 
the two-sphere results of Reed and Morrison 
(10) to estimate the electrophoretic velocity 
of a test sphere of one type in a dilute suspen- 
sion of spheres of different zeta potentials. 
These studies indicate that particle-particle 
interactions are relatively weak; detailed ex- 
perimental work (12) shows that this is true 
for homogeneous systems. 

In our paper on homogeneous suspensions 
(12), mobilities of individual tracer particles 
were measured over a wide range of particle 
concentrations and ionic strengths. This in- 
volved the use of red blood cell ghosts and 
tracer particles with almost identical electro- 
kinetic properties. With these suspensions, the 
electrical conductivity decreased rapidly with 
increasing cell concentration, indicating strong 
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electrostatic interactions between particles. 
Nevertheless, the mobilities of individual par- 
ticles were relatively insensitive to the presence 
of neighbors. The speed of translation varied 
as ( 1 - C), C being the volume fraction of 
particles, up to the highest volume fractions 
achieved. Since this 0(C ) term arises solely 
from the backflow necessary to conserve vol- 
ume in a space filling suspension, particle- 
particle interactions cancel insofar as electro- 
phoresis is concerned. The particles behave as 
electromechanical, as well as optical, ghosts. 
The purpose of the current paper is to report 
on the behavior of heterogeneous suspensions, 
i.e., those containing particles of differing mo- 
bilities. 

For this study we needed suspensions of 
stable, transparent, colloidal particles with dif- 
ferent electrophoretic mobilities, along with 
tracer particles with the same electrokinetic 
properties. These can be prepared from sus- 
pensions of human red blood cells (RBCs) via 
osmotic lysis (13), which yields transparent 
spherical ghosts with electrophoretic mobilities 
equal to those of the native cells (12). The 
ghosts and native RBCs were fixed with glu- 
taraldehyde to render them “rigid” and stable 
to reduction in ionic strength (14). Treatment 
of fixed cells with Neuraminidase (15) reduces 
the electrophoretic mobilities of RBCs and 
ghosts substantially below those of native cells. 
Heterogeneous suspensions were prepared by 
mixing untreated ( or fast ) ghosts with enzyme 
treated (or slow) ghosts. 

Our experiments show that the electropho- 
retic mobilities of either fast or slow tracers 
are relatively insensitive to the presence of 
other particles in the suspension. Tracer mo- 
bilities were linear functions of total suspen- 
sion volume fraction out to the point where a 
change in rheological behavior occurs. At this 
volume fraction (approximately 0.45 at the 
ionic strength used here), the particles ap- 
peared to “lock together” and relative move- 
ment of fast and slow particles was inhibited. 

The remainder of this paper is divided into 
four sections. Section II deals with the exper- 
imental details and Section III with the results. 


Section IV discusses previous electrokinetic 
experiments and existing models. Concluding 
remarks are in Section V. 

II. EXPERIMENTAL TECHNIQUES 

The ghost preparation procedures used in 
this work closely follow the protocol laid out 
by Zukoski and Saville (12). Fresh human 
blood was drawn by venipuncture into 5-ml 
Vacutainers ( Vacutainer Systems, Rutherford, 
NJ ) and anticoagulated with 2. 1 mg disodium 
EDTA per milliliter of blood. 

Ghosts were prepared by an iso-ionic lysis 
technique described by Billah et al. (13). 
Whole blood was washed three times in 0. 1 5 
M PBS (containing 2.54 X 10 -2 MNa 2 HP0 4 , 
4. 62 X 1 0 ' 3 M NaH 2 P0 4 , and 0. 1 20 M NaCl, 
pH 7.5 ± 0.2) and one time in 0.15 M PBS 
+ MgCl 2 (the same 0.15 M PBS buffer de- 
scribed above but containing 2 X 10" 3 Mf 
liter MgCl 2 ). Packed cells were resuspended 
to a volume fraction of 0.25 with 0. 1 5 M PBS 
+ Gly (0. 1 5 M PBS + MgCl 2 containing 2 M / 
liter ethylene glycol). The suspensions were 
incubated 15 min at room temperature and 
then 15 min at 0°C. Twenty parts of ice cold 
0. 1 5 M PBS + MgCl 2 was added rapidly to the 
chilled suspension for every part of original 
cell pack volume (final cell volume fraction, 
0.04 ) . Lysis was evident by clearing of the tur- 
bid appearance. The ghost cells were packed 
by centrifugation at 1 6,000 g for 20 min in a 
Beckman L5-65 Preparative Ultra-Centrifuge 
using a type 2 1 rotor. The resulting cell pack 
was washed once in 0. 1 5 M PBS + MgCl 2 and 
the lysis procedure repeated until the desired 
amount of hemoglobin was removed from the 
ghosts; four lysis steps were usually required 
to produce hemoglobin free ghosts. 

The ghosts are uniformly spherical and, as 
reported by Billah et al. ( 1 3 ) , impermeable to 
large macromolecules such as dextrans or al- 
bumin. Little evidence for ghost fragmentation 
was observed. As reported earlier (12), the av- 
erage ghost volume was 120 ± 20 j/m 3 , a °d 
neither average size nor the size distribution 
was altered by fixation. 
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All electrokinetic studies reported in this 
paper were carried out with fixed ghosts and 
fixed RBCs. Fixation stabilized the suspen- 
sions against degradation over time, as well as 
the action of the enzyme treatment described 
below, and rendered the cells rigid. Fixation 
was carried out by washing packed RBCs or 
ghosts twice in 0. 15 M PBS and resuspending 
them to a final volume fraction less than 0.10 
in 0.15 M PBS containing 1% (w/v) glutar- 
aldehyde (Sigma, electron microscopy grade). 
The RBC suspensions were checked for he- 
moglobin release and discarded if any was 
present. This fixation process, adapted from 
the work of Vassar et al. (14), produced sta- 
bilized cells that retained their original shape 
and size distribution. The fixed RBCs showed 
no signs of rupture when suspended in deion- 
ized water and viewed under the microscope. 
The fixation process did not alter the size or 
shape distribution of the ghost cells, and the 
ghost cells were stable when washed into media 
of lower ionic strength. However, if the process 
was reversed, the cells appeared to collapse 
under the action of the external osmotic pres- 
sure. Consequently, fixed ghosts were always 
washed into lower ionic strength media. 

Neuraminidase treatment followed the 
work of Nordt et al. { 15 ). Type V neuramin- 
idase (Sigma) was reported by the manufac- 
turer to have a specific activity of 4.9 units/ 
mg protein using iV-acetyl neuraminic-lactose 
as a substrate and 0.5 units /mg protein with 
bovine submaxillary mucin as a substrate. One 
unit will liberate 1 .0 /^mole of iV-acetyl neur- 
aminic acid per minute at pH >5.0 at 37°C us- 
ing the indicated substrate. Fixed RBCs and 
ghosts were washed three times in 0. 1 5 M PBS 
and finally resuspended to a volume fraction 
of approximately 0.20. Neuraminidase was 
dissolved in 0.15 M PBS and added to sus- 
pensions such that the final concentration was 
0.27 units per 10 10 cells. Then the suspensions 
were incubated for 2.5 h at 37 °C, washed two 
times in 0.5 M PBS, and stored at 4°C until 
use or resuspended in 0. 1 5 M PBS containing 
1% (w/v) glutaraldehyde. The shape and size 
distributions of the fixed RBCs and fixed 


ghosts were not altered by this enzymatic 
treatment. The enzyme treated cells showed 
signs of aggregation in 0. 15 M PBS but were 
easily redispersed under mild shear. Before and 
after treatment the mobilities of RBCs and 
ghosts were the same to within experimental 
error (5%). Enzyme treatment reduced the 
RBC and ghost mobilities by almost 70% (Ta- 
ble I ) . Cells were treated in amounts such that 
all electrokinetic experiments reported here 
could be carried out using RBCs and ghosts 
from the same neuraminidase preparation. 
Experiments carried out using another enzyme 
preparation, however, showed similar results. 

Electrokinetic studies were carried out using 
a Rank Mark II Electrophoresis Apparatus 
(Rank Brothers, Bottisham, England) 
equipped with silver/ silver chloride electrodes 
and a flat cell. The silver/ silver chloride elec- 
trodes eliminated electrode polarization, 
which would occur with platinum or palla- 
dium electrodes at the ionic strengths of these 
experiments. Conductivity measurements 
were made on suspensions in a flat quartz 
electrophoresis chamber with an Altex (Beck- 
man ) conductivity bridge capable of measur- 
ing conductivities at 80 and 1000 Hz. Within 
this frequency range the conductivity of sus- 
pensions did not vary. All measurements were 
carried out at 25 ± 0.02°C. 

To measure mobilities, small amounts of 
RBC tracers were added to suspensions of 
ghosts to give a RBC tracer concentration near 
that used in analytical electrophoresis. Then 
the electrophoresis chamber was filled with the 
suspension just prior to making the mobility 
and conductivity measurements. After the 
measurement, small volumes of the suspen- 
sion (0.5- 1.0 ml) were removed for subse- 
quent volume fraction determination and re- 
placed with fresh buffer. This process was re- 
peated until the ghost volume fraction reached 
a low value, typically less than 0.08. The 
tracer/ghost volume fraction ratio stayed ap- 
proximately constant during the dilution pro- 
cess. 

As discussed earlier (12), the velocity of 
particles in a dilute system is a sum of the 
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particle’s electrophoretic velocity and the 
electroosmotic velocity of the fluid in the 
chamber. The electrophoretic velocity of a 
particle is relatively independent of position, 
whereas the electroosmotically driven flow is 
not. Moreover, if the volume fraction is low, 
a parabolic flow profile exists. In such situa- 
tions the mobility can be measured at the sta- 
tionary plane where the electroosmotic veloc- 
ity of the fluid is zero. This position depends 
on chamber geometry and can be found a 
priori. On the other hand, if the suspension is 
concentrated, the position of the stationary 
plane will depend on the constitutive relation- 
ship between stress and strain, as well as 
chamber geometry. Under these circum- 
stances, the stationary plane cannot be found 
in an a priori fashion. Nevertheless, the elec- 
trophoretic mobility can be separated from the 
effects of bulk flow since the net flow through 
the chamber is still zero. Integrating the tracer 
velocity profile across the chamber yields the 
electrophoretic velocity directly. As discussed 
by Hunter (16), this technique reduces errors 
associated with the position of the stationary 
level even under standard analytical electro- 
phoresis conditions. 

In the experiments described here, the ve- 
locity of individual tracer particles was mea- 
sured at 10-15 positions across the chamber, 
starting close to one wall and making the final 
measurements close to the opposite wall. 
However, to avoid excluded-volume effects, 
positions very close to either wall were 
avoided. At each depth, between 4 and 10 
tracer velocities were measured in alternately 
opposite directions. The resultant velocity 
profile was integrated numerically to calculate 
the average tracer velocity and tracer mobility. 
Typical field strengths were less than 1 0 V / 
cm, but at the elevated vokime fractions where 
the suspensions were highly viscous, field 
strengths of 20 V/cm were used. Nevertheless, 
the average tracer velocity always varied lin- 
early with the applied field strength. 

The fixed ghost volume fractions were de- 
termined from the suspension number con- 
centration and the average volume of the ghost 


particles. Spherical ghosts prepared by the 
ethylene glycol lysis technique (13) could be 
counted and sized accurately using a Particle 
Data Elzone Celloscope (Particle Data, Elm- 
herst, IL) calibrated with standard latex par- 
ticles. This electronically determined size dis- 
tribution was checked against a microscopi- 
cally determined size distribution (12). To 
within experimental error, the two gave the 
same distribution and average size. 

In experiments where a 50-50 mixture of 
fast and slow ghosts was used, dilute suspen- 
sions of each were made up and adjusted to 
the same number concentration. Equal vol- 
umes of each suspension were mixed, and the 
resulting suspension packed by centrifugation. 
The difference between the concentrations of 
fast and slow ghosts was less than 1%. 

Most of our electrokinetic studies were car- 
ried out in 0.045 M PBS (0.15 M PBS diluted 
with deionized water). Here the suspensions 
were washed two times in 0.045 M PBS and 
allowed to stand overnight at 4°C in 0.045 M 
PBS prior to a final washing. This assured 
equilibration of ionic strength across the 
membrane. 

III. RESULTS 

Some preliminary experiments were carried 
out on homogeneous suspensions of native 
tracers and ghosts and enzyme treated tracers 
and ghosts in 0.15 M PBS. At this ionic 
strength, the enzyme treated ghosts ( and trac- 
ers) have electrophoretic mobilities of —0.56 
(n m/s)/(V/cm), which corresponds to a 
surface potential of approximately —7 .0 mV . 
These suspensions were very viscous and non- 
Newtonian. This suggests that the enzyme 
treated cells were colloidally unstable under 
these conditions and coagulation caused the 
anomalous results. All the results reported here 
are for systems where the electrolyte was 0.04f 
M PBS and the particles were stable. 

Figures 1 and 2 show conductivities for ho- 
mogeneous suspensions in 0.045 M PBS 
Conductivities for homogeneous suspensions 
of fast and slow ghosts coincide, to within ex- 
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Fig. 1. The conductivity of suspensions of fixed un- 
treated ghosts as a function of suspension volume fraction. 
The background electrolyte is 0.045 M PBS (pH 7.4). 

perimental uncertainty at volume fractions as 
large as 0.6, suggesting that enzyme treatment 
does not significantly alter membrane prop- 
erties at this ionic strength. Due to the high 
ionic strength and low surface charge in these 
suspensions, changes in the surface charge 
have a negligible effect on the bulk conductiv- 
ity. The deviation of the limiting slope from 
the value of — 1.5 given by Maxwell’s equation 
may be due to membrane or double-layer 
conductance. Figure 3 shows the electropho- 
retic mobility of fast tracers in suspensions of 
fast ghosts, along with the mobility of slow 
tracers in slow ghosts. The tracer mobility var- 
ies linearly with ghost volume fraction, inde- 
pendent of enzyme treatment. For these sus- 
pensions, the electrophoretic increments de- 
rived by linear regression are —0.97 for 
suspensions of fast tracers in fast ghosts and 
- 1 .00 for suspensions of slow ghosts. 3 

Figure 4 depicts the mobility-volume frac- 
tion relationships for fast tracers in suspen- 
sions of slow ghosts and slow tracers in sus- 
pensions of fast ghosts. There is a linear vol- 
ume fraction dependence of fast tracer 
electrophoretic mobilities in slow ghost sus- 
pensions, and the electrophoretic increment 
is, to within experimental error, the same as 
the electrophoretic increment for fast tracers 
suspended in fast ghosts. Similarly, the mo- 

3 The electrophoretic increment, A/x, is defined through 
an expansion of the particle’s mobility, n e , in total 
suspension volume fraction C, as n e = M«(l + AmC 
4- • • • ), where is the mobility as C -*■ 0. 



VOLUME FRACTION(C) 


Fig. 2. The conductivity of suspensions of fixed enzyme 
treated ghosts as a function of suspension volume fraction. 
The background electrolyte is 0.045 M PBS (pH 7.4). 

bility of slow tracer particles in suspensions of 
fast ghosts is linearly dependent on ghost vol- 
ume fraction out to a volume fraction of 0.4. 
The electrophoretic increment for the linear 
portion of the curve is —1.15; within experi- 
mental uncertainty, this equals the electro- 
phoretic increment for slow tracers in slow 
ghost suspensions. 

Above a volume fraction of 0.4 the suspen- 
sions become substantially non-Newtonian 
and the electrophoretic mobility of the slow 
tracers in the fast ghosts increases. At the 
highest volume fraction studied (0.59), the 
velocity profile had a blunt shape at the cham- 



VOLUME FRACTION (C) 


FIG. 3. The electrophoretic mobility of fixed untreated 
RBC tracers in suspensions of fixed untreated ghosts (O) 
and the electrophoretic mobility of fixed enzyme treated 
tracers in suspensions of fixed enzyme treated ghosts (□) 
as a function of suspension volume fraction. The back- 
ground electrolyte is 0.045 M PBS (pH 7.4). 
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Fig. 4. The electrophoretic mobility of fixed untreated 
RBC (fast) tracers in suspensions of fixed enzyme treated 
(slow) ghosts (O) and the electrophoretic mobility of fixed 
enzyme treated (slow) tracers in suspensions of fixed un- 
treated ghosts (□) as a function of suspension volume 
fraction. The dashed line represents a mobility that varies 
as ( 1 - C). The background electrolyte is 0.045 M PBS 
(pH 7.4). 

ber center and the suspension displayed some 
elasticity. Here the slow tracer mobility 
equaled that expected for fast tracers in fast 
ghost suspensions. This suggests that the slow 
tracer particles were locked in place and car- 
ried along with the fast ghosts. It should also 
be noted that the mobility intercept for fast 
tracers shown in Fig. 4 is lower than those re- 
ported in Figs. 3 or 5. For this batch of fast 
tracers, the mobilities decreased as a result of 
fixation and, consequently, these tracers were 
used only in the experiments shown in Fig. 4. 

Figure 5 shows electrophoretic mobilities for 
fast and slow tracer particles plotted as a func- 
tion of volume fraction for a 50-50 mixture 
of fast and slow ghosts. Up to volume fractions 
of 0.35, the mobilities of both fast and slow 
tracer particles are linearly dependent on total 
suspension volume fraction. Moreover, the 
electrophoretic increments for both the slow 
and the fast tracer particles are, to within ex- 
perimental error, the same as for fast tracers 
suspended in fast ghosts and slow tracers sus- 
pended in slow ghosts. 

At volume fractions greater than 0.36, a 
large amount of time dependent structure was 


observed to build up in the heterogeneous sus- 
pensions. This structure took the form of wavy 
alternating dark and light bands perpendicular 
to the direction of particle motion and moving 
with the particles. At a volume fraction of 0.40, 
the structure was so intense as to make mea- 
surements irreproducible. Here the tracer par- 
ticles gradually slowed down at the center of 
the chamber until the entire system moved in 
plug flow. Reversing the field caused the struc- 
ture to break down over a period of 8 to 10s 
until the tracers appeared suspended in a uni- 
formly grainy background. Suspension struc- 
ture reformed moving in the opposite direction 
within approximately 30 s. During this time, 
the tracer particles would slow down and form 
a flat velocity profile. In the volume fraction 
range from 0.36 to 0.30, structure became less 
evident, and tracer velocities were independent 
of the length of time the electric field had been 
applied. No structure was observed at volume 
fractions less than 0.30. It should be noted that 
with homogeneous suspensions, little or no 
structure was observed at volume fractions as 
large as 0.59. At this volume fraction, a small 
amount of elasticity was observed, but after a 
short time (10-15$) the tracer velocities were 
independent of time. The data presented in 



Fig. 5. The mobility of fixed untreated RBC (fast) trac- 
ers (O) and fixed enzyme treated tracers (□) in suspensions 
of a 50-50 mixture of fixed untreated and fixed enzyme 
treated ghosts. The dashed line represents a mobility that 
varies as ( 1 - C). The background electrolyte is 0.045 M 
PBS (pH 7.4). 
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Fig. 5 were taken on suspensions only where 
steady mobilities were observed. The inability 
to gather meaningful tracer velocity profiles 
above volume fractions of 0.4 in 50-50 mix- 
tures of slow and fast ghosts suggests that a 
mechanism similar to the “locking-up” ob- 
served for slow tracers in suspensions of fast 
ghosts was present. 

The electrophoretic data presented in Figs. 
3-5 are rather surprising. If the suspension 
volume fraction is kept below a certain value, 
the electrophoretic mobilities of the tracer 
particles are almost independent of the elec- 
trokinetic properties of the other particles 
which make up the bulk of the suspension 
volume fraction and depend only on the over- 
all particle concentration. 

III. DISCUSSION 

Previous studies have established sialic acid 
as the main charge determining moiety on the 
red blood surface (17) and treatment of RBCs 
with neuraminidase is known to release vir- 
tually all external membrane-bound sialic acid 
(15). In 0.15 M NaCl (pH 7.2), the electro- 
phoretic mobility of fresh neuraminidase 


treated RBCs is virtually zero, suggesting a 
complete removal of the red cell’s electroki- 
netic charge. Subsequent glutaraldehyde fix- 
ation of RBCs that had been treated with 
neuraminidase to reduce the mobility almost 
to zero increases cell electrophoretic mobilities 
to -0.45 ± 0.09 (^m/s)/(V/cm) (14). The 
source of this increase in mobility is not 
known; however, the increase is stable over 
time. Fixation of native RBCs and native RBC 
ghosts with glutaraldehyde raises the cell mo- 
bility from —1.08 ± 0.05 (yum/s)/(V/cm) to 
-1.26 ± 0.07 ( # tm/s)/(V/cm) in 0.15 M 
NaCl, pH 7.2. Again this increase is of un- 
known origin and is stable with time (15). 
Table I gives the electrophoretic mobilities at 
infinite dilution of the native and enzyme 
treated particles used in this study. According 
to the table, neuraminidase treatment reduced 
the mobilities of the RBCs and ghosts by al- 
most 70%. 

For the studies at hand, particles of sub- 
stantially different mobility were required, and 
the enzyme treatment described here provided 
these particles. At the ionic strengths used here, 
fixed RBCs and fixed ghosts had the same 
electrophoretic mobilities as did the enzyme 


TABLE I 


Electrokinetic Properties 


Sample* 

Aa 6 

a c 

((im/sKV/cm) 

r(mV) 

Untreated RBC/ghosts in 
0.15 M PBS 

-1.26 

0.1 L0_ . 

-1.26 ±0.06 

-16.2 

Enzyme-treated RBC/ghosts • 
in 0.15 M PBS 

-1.90 

_ 

-0.56 ± 0.08 

-7.0 

Untreated RBC/ghosts in 
0.045 M PBS 

-1.14 

0.170 

-2.20 ± 0.08 

-29.0 

Enzyme-treated RBC/ghosts 
in 0.045 M PBS 

-1.17 

0.170 

-0.80 ±0.06 ' 

-10.5 


a Note that all the cells were fixed with glutaraldehyde. 

b The conductivity increment, Act, was calculated from the limiting slope of the conductivity-volume fraction re- 
lationship for suspensions of ghosts alone. 

f The relative conductivity of the particles in suspension, a, is calculated from Act and is valid only for the ghost 
particles and does not necessarily hold for the tracer RBC ( 12). 

d The electrophoretic mobility, measured at infinite particle dilution, was the same for RBC tracers and ghosts 
for both untreated and enzyme treated samples. 

f The electrokinetic surface potential, f, as calculated from the Smoluchowski formula. 
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treated fixed RBCs and enzyme treated fixed 
ghosts. Moreover, the enzyme treatment did 
not substantially alter the distribution of 
electrophoretic mobilities, which remained 
roughly Gaussian with a standard deviation 
of ±5% (17). Thus, enzyme treated RBCs 
faithfully duplicated the behavior of the en- 
zyme treated ghost particles. 

The electrokinetic behavior of the homo- 
geneous suspensions used in this work was 
similar to that found in the earlier study (12). 
Figure 3 shows that the mobilities are linear 
in the volume fraction of particles and the 
electrophoretic increment is, to within exper- 
imental error, — 1.0. This sort of behavior can 
be attributed to the irrotational nature of the 
velocity field of identical particles undergoing 
electrophoresis when the diffuse layer is thin, 
as it is here. The argument hinges on the fact 
that the electrostatic and hydrodynamic in- 
teractions between a pair of particles cancel 
and goes as follows. 

The velocity of a pair of identical spheres 
undergoing electrophoresis can be calculated 
using the method of reflections. Starting with 
the velocity of an isolated sphere, the velocity 
of the other sphere is calculated, ignoring its 
effect on the velocity of the first. Then the ve- 
locity of the first sphere is corrected for the 
presence of the second, after which the velocity 
of the second is recalculated using the cor- 
rected velocity of the first sphere. When re- 
peated ad infinitum, this yields the velocity of 
either sphere in the presence of the other. 

First we note that a single particle translating 
through an infinite fluid produces a certain 
velocity field, u(r), in the surrounding fluid. 
A second particle placed in this fluid will re- 
spond to the viscous and electrical forces. The 
velocity of the second particle differs from 
what it would be if it were all alone due to the 
shear flow and alteration in the electric field 
produced by the first particle. The velocity of 
a single particle placed in a shear flow can be 
found from Faxen ’s law ( 1 8 ) as 

U = U 0 + u(r) + V 2 u(r). [1] 

Here Uo represents the electrophoretic velocity 


of the second sphere, i.e., the product of the 
mobility and the local field at the particle’s 
location, viz., 

U 0 = Me(E + SE). [2] 

E is the undisturbed field, 5E is the disturbance 
to the uniform field, and ^ is the electropho- 
retic mobility of the isolated particle. Note that 
Faxen ’s law does not account for the effect of 
the second particle on the velocity of the first 
and, at this stage, Eq. [1] simply represents the 
uncorrected velocity of the second particle in 
the presence of the first. To complete the first 
stage of the calculation we need the velocity 
field and the deformation to the electric field 
produced by the first particle. For a single par- 
ticle undergoing electrophoresis, it is easily 
shown that when the diffuse layer is thin the 
velocity field around it is irrotational and given 
by 

u(r) = -Me^- 3 [I - 3r~ 2 rr] • E. [3] 

Here r is the position vector scaled on the par- 
ticle radius and I is the unit tensor. The dis- 
turbance to the electric field produced by the 
first sphere is 

5E(r) = ^r -3 [ I — 3r~ 2 rr] • E. [4] 

Combining Eqs. [l]-[4] shows that in the first 
approximation the second particle moves at 
the speed it would if the first particle were ab- 
sent. At the next stage in the process we cal- 
culate the correction to the velocity of the first 
particle due to the presence of the second; it 
too vanishes. Continuing the process shows 
that the two identical spheres move as if they 
were all alone, regardless of their orientation 
with respect to the field. This is the result found 
by Reed and Morrison (10) by means of a 
more detailed calculation. To calculate the 
velocity of a test particle in a space filling sus- 
pension requires an accounting for: (i) the 
backflow due to the motion of all other par- 
ticles and (ii) the alteration in the strength of 
the field due to the particles. The former cal- 
culation is straightforward for systems where 
the particles are distributed at random around 
the test particles (18) and yields 
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Me = Me( 1 ” C) [5] 

for the mobility. It may also be noted that the 
argument given above holds for concentrated 
suspensions where, due to the irrotational na- 
ture of the velocity and electric field, the hy- 
drodynamic and electrostatic effects cancel up 
to particle concentrations where the diffuse 
layers overlap. 

The theory of electrophoresis in a dilute, 
homogeneous system was worked out by An- 
derson ( 1 1 ) by adapting Batchelor’s results for 
sedimentation (18). However, the fact that the 
structure of the mean field is altered by the 
particles in a space filling suspension is absent 
from Eq. [ 5 ] and, according to Anderson (19), 
accounting for this effect with nonconducting 
particles adds —j^C to the equation. Nev- 
ertheless, the experimental results clearly show 
that the electrophoretic increment for homo- 
geneous suspensions is close to —1, not 
although the reason for this discrepancy is un- 
known. It could be due to the slightly poly- 
disperse, slightly nonspherical nature of the 
particles. Or there might be a subtle error in 
the theory. In any event, the problem deserves 
further study. 

The results from experiments where the 
mobility of the tracer differs from that of the 
suspension particles were shown in Fig. 4. Here 
the electrophoretic increments derived from 
linear regression of the experimental data are 
- 1 .03 and — 1. 1 5 for fast tracers in slow ghosts 
and slow tracers in fast ghosts, respectively. 
These are not significantly different from —1, 
as the dashed line in the figure clearly indicates. 
With fast or slow tracers in slow or fast ghosts, 
the distribution of neighbors around the tracer 
is nonuniform and the increment should differ 
from — 1 . Evidently the effect of the nonuni- 
formity is weak. 

The situation with heterogeneous suspen- 
sions shown in Fig. 5 is even more interesting. 
Here the increments are, to within the accu- 
racy of the measurement, — 1 for either sort of 
tracer in 50-50 mixtures of fast and slow 
ghosts. To compare this result with a theory 
of the sort embodied in Eq. [5] requires 


knowledge of particle distribution around the 
test particle. 4 The methodology has been 
worked out for sedimentation in dilute systems 
( 20-22 ) and adapted for electrophoresis ( 23 ) . 
For binary systems the result is 

A*e,/ = Ml[l -C-0.13^Q] [6] 

to O(C) . Here the i, j subscripts identify the 
particle type and 

Vjj = 1 - 17 ] 

The numerical coefficient 0. 1 3 depends on the 
pair distribution, and its value is extremely 
sensitive to the details of the particle-particle 
interaction at close separations. Thus, the nu- 
merical value reported here may be revised 
when more accurate calculations are available. 
Upon using the results from Table I we have 

- 1-04C] [8] 

and 

Me.s = Me,s[I ~ 0.89C] [9] 

for fast and slow particles in 50-50 mixtures 
of the two. The theory indicates, as expected, 
that fast particles are slowed somewhat by in- 
teractions with slower neighbors (Eq. [8]). 
Similarly, slow particles are helped along by 
interactions with faster neighbors (Eq. [9]) 
compared with their behavior in homogeneous 
systems. Nevertheless, the effect is weak and 
not very different from the experimental re- 
sults in Fig. 5. Finally, it should be noted that 
these calculations are for dilute systems, and 
multiparticle interactions have been ignored. 

V. CONCLUSIONS 

The experimental work reported here on 
heterogeneous systems shows that particle- 
particle interactions in suspensions are sur- 
prisingly weak insofar as electrophoretic pro- 
cesses are concerned. With homogeneous sys- 
tems it can be shown that electrostatic and 
hydrodynamic interactions cancel so the ab- 

4 Anderson’s calculation for heterogeneous systems (II) 
is applicable only when the spheres around the test particle 
are “well stirred.” 
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sence of particle-particle interactions is readily 
understood, even in concentrated systems. 
However, with heterogeneous systems inter- 
actions do not cancel, and the distribution of 
particles around a test particle is no longer 
uniform in a dilute system. However, although 
stronger interactions might have been antici- 
pated, the redistribution of particles has a 
rather weak effect on electrophoretic mobility. 
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Two problems of similar mathematical structure are studied : the thermocapillary 
motion of bubbles and the electrophoresis of colloidal particles. The thermocapillary 
motion induced in a cloud of bubbles by a uniform temperature gradient is 
investigated under the assumptions that the bubbles are all the same size, that the 
surface tension is high enough to keep the bubbles spherical, and that the bubbles are 
non-conducting. In the electrophoresis problem, the particles, identical spheres 
having a uniform zeta potential, are suspended in an electrolyte under conditions 
that make the diffuse charge cloud around each particle small when compared with 
the particle radius. For both problems, it is shown that in a cloud of n particles 
surrounded by an infinite expanse of fluid, the velocity of each sphere under creeping 
flow conditions is equal to the velocity of an isolated particle, unchanged by 
interactions between the particles. However, when the cloud fills a container, 
conservation of mass shows that this result cannot continue to hold, and the average 
translational velocity must be calculated subject to a constraint on the mass flux. 
The computation requires ‘renormalization but it is shown that the renormalization 
procedure is ambiguous in both problems. An extension of Jeffrey’s (1974) second 
group expansion, together with the constraint of conservation of mass, removes the 
ambiguity. Finally, it is shown that the average thermocapillarv or electrophoretic 
translational velocity of a particle in the cloud is related to the effective conductivity 
of the cloud over the whole range of particle volume fractions, provided that the 
particles are identical, non-conducting and, for the thermocapillarv problem, 
inviscid. 


1. Introduction 

If a temperature gradient is applied to a viscous fluid containing bubbles, the 
bubbles move towards the hotter fluid, owing to the dependence of surface tension 
on temperature. Similarly, if an electric field is applied to an electrolytic solution 
containing charged particles, the particles move in directions dictated by their 
charge. Under conditions described below, the problems of thermocapillarv motion 
and electrophoresis are mathematically similar, the governing equations being the 
same except for one boundary condition. We wish to calculate the average velocity 
of a particle in a cloud (we shall understand the word ‘particle ’ to encompass bubbles 
as well) when the cloud is placed in a uniform gradient of temperature or electric 
potential and the particle concentration is low. The problem has several features in 
common with that of calculating the average settling velocity of a particle in a 
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sedimenting cloud, a problem that was the subject of a paper by Professor Batchelor 
(Batchelor 1972) that gave a strong stimulus to research on the mechanics of 
suspensions. Its significance lay not just in its calculation of a sedimentation 
velocity, but equally in the way it addressed wider themes, such as the use of 
statistical methods to calculate average properties of suspensions and the solution of 
a subtle technical problem, namely how to express the average properties of a 
suspension of particles in terms of convergent integrals. The success of this and later 
work on the mechanics of suspensions, together with parallel progress on other flow 
systems containing small lengthscales, lead Batchelor (1976) to identify a new 
branch of specialization in fluid mechanics: microhydrodynamics. 

One of the technical difficulties that must be solved during a calculation of the bulk 
properties of a suspension is the appearance of non-convergent integrals in various 
formal expressions for the quantities of interest, such as the average sedimentation 
velocity of the particles or the effective viscosity of a suspension. Earlier workers had 
already noted and tackled this difficulty (Burgers 1942; Pyun & Fixman 1964; see 
also references in Jeffrey 197 7 ) and so its widespread occurrence was known when 
Batchelor (1972) published his method for overcoming it. Because of this, the new 
method, later called renormalization, was quickly applied to other problems (for 
example, Batchelor & Green 1972; Jeffrey 1973; Chen & Acrivos 1978). In addition, 
the method attracted a great deal of examination for its own sake, and was later 
refined and extended in many different ways (Jeffrey 1974; Willis & Acton 1976; 
McCoy & Beran 1976; Goddard 1977 ; O’Brien 1979; Felderhof, Ford & Cohen 1982). 

Several concerns can be identified in the papers that reanalysed Batchelor’s 
method of renormalization. The first was whether renormalization was of physical or 
mathematical significance. It must be remembered that both Rayleigh (1892) and 
Einstein (1906) had evaluated non-convergent sums or integrals in the course of 
deriving firmly established results, and consequently non-convergence seemed to 
some to be a purely mathematical issue and to others a detail requiring just some 
additional physical justification. Thus the first notable feature of the new approach 
was its stress on the necessity of renormalization. In this regard, an important benefit 
of tackling the sedimentation problem first (rather than, say, the viscosity problem) 
was the fact that the integrals were always infinite, and a discussion of ‘physically 
significant ’ evaluations was impossible. 

A second concern, closely related to the first, was a desire to simplify a key step 
in the original analysis which entailed looking ‘for a quantity whose mean is known 
exactly from some overall condition or constraint in the specification of the problem 
and whose value at at 0 [the position of a test particle - notation similar to Batchelor’s 
is defined in §3] has the same long-range dependence on the presence of a sphere at 
jt 0 + r as the velocity of the test sphere; and, once found, the difference between V 
[the sedimentation velocity] and the mean of this quantity can be expressed as an 
integral... which can be legitimately... evaluated explicitly. ’ This description and the 
finding of renormalizing quantities were not easy to grasp, and even rereading 
Batchelor (1972) today, one is struck by the difficulty of some of the ideas; who, for 
example, would find it obvious to choose the divergence of the deviatoric stress for 
use in Batchelor’s equation (3.10)? 

In fact, a discussion of the physical differences between a cloud of particles that 
fills a container and one that is surrounded by clear fluid appears as early as the 
second paragraph of Batchelor (1972), and this is the essence of a physical 
interpretation. Nevertheless, renormalization was presented as a means of removing 
a convergence problem which had a physical interpretation, rather than as a means 
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of accounting for backflow that had a mathematical formulation. Similarly, although 
it was clearly expected that the renormalizing quantity would be physically related 
to the integrand being renormalized, it was not until Jeffrey (1974) that a 
methodology was introduced that enforced this requirement systematically. The 
inadequacy of the purely mathematical interpretation was clearly displayed when a 
problem was discovered in which several choices of the renormalizing quantity 
appeared to be possible (Chen & Acrivos 1978), a problem we discuss below. 

Another reason why Batchelor s method w r as re-examined stemmed from issues 
connected with the distinction between the value of the external field applied to the 
suspension and the mean value of that field within the suspension. Ever since the 
investigation of effective viscosity b}’ Einstein (1906), the possible difference between 
the applied field — in other words, the field ‘at infinity ’ - and the mean field had been 
a point of debate. Investigators studying the electrostatic and elastostatic properties 
of dispersions of particles could formulate problems in which the suspension was 
surrounded by clear fluid or matrix, and thereby they could separate the applied and 
the mean fields. For sedimentation, however, the very nature of the problem made 
it necessary to deal with a space-filling suspension from the outset, leaving these 
issues unresolved. They were addressed in subsequent work (cited above) and the 
results obtained by Batchelor verified ; and although the original method has now 
been filled out in completeness and rigour, it has yet to be replaced by an approach 
that is simpler or more direct, at any rate for problems to which it is suited. 

Batchelor's 1972 paper stimulated work in other ways as well. It drew attention 
to the difference between the sedimentation of a random dispersion of freely moving 
spheres, the sedimentation of a periodic array of spheres and the flow through a 
random array of fixed spheres. The last case is a ‘strong interaction’ problem that 
cannot be solved by Batchelor's method, and the efforts to understand this and find 
an alternative method lead to the averaged-equation approach (Childress 1972; 
Saffman 1973; Howells 1974; Hinch 1977). It is interesting to note that the new 
approach had first to explain Batchelor’s method in its own terms (Hinch 1977). 
Finally, it is instructive to list the topics not treated in the original paper which 
received attention in later years. For example, the paper contains comments on the 
importance of particle motions for the probability distribution for the particle 
positions, even though in the actual calculation, specific distributions were simply 
assumed. This lead was taken up by many later papers, and the motions of pairs of 
particles under the influence of hydrodynamic and non-iiydrodynamic forces have 
been explored in detail. Other extensions to the original investigation include the 
addition of non-hydrodynamic effects, such as Brownian motion, the effects of 
unequal sizes and of inhomogeneous conditions. 

It is a tribute to Batchelor's idea that a new problem should be discovered that 
leads us to re-examine the basis of his method. This problem, which concerns the 
thermocapillary motion or electrophoresis of particles, brings a new element into the 
analysis and adds another example of a well-known difficulty. The new element 
concerns the ‘applied field’. Specifically, when Jeffrey (1974) presented a general 
treatment of Batchelor's renormalization, all the problems known at that time could 
be summarized symbolically in terms of a quantity S to be averaged and a quantity 
/f-the applied field - whose average was known. In the thermocapillary and 
electrophoretic problems, however, there are two quantities whose averages are 
specified and which determine the velocities of the particles, because both the applied 
field (electric or thermal) and the volume flux across any boundary spanning the 
suspension are known and can be used as renormalizing quantities. The well-known 
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difficulty is the possibility, first discovered by Chen & Acrivos (1978). that a studv 
of two-particle interactions may not be sufficient to prove which of several choices 
gives the correct renormalization strategy. 

The study of thermocapiilarv phenomena has taken on new significance with the 
de\ elopment of the space shuttle and the opportunities for experimenting and 
manufacturing under near-weightless conditions. For example, unwanted gas 
bubbles in a mixture are usually removed by buoyancy forces, but this force is not 
available in an orbiting laboratory. In fact, even on Earth a viscous fluid can retain 
bubbles for a significant length of time if the bubbles are small enough, and under 
near-weightless conditions forces other than buoyancy must be found to move 
bubbles, which is why thermocapillary forces are of interest. The motion of a single 
bubble was studied by Young, Goldstein & Block (1959) and the motion of two 
bubbles b\ Mev\appan, \\ ilcox & Subramanian (1983), Meyvappan & Subramanian 
(1984), Anderson (1985a) and Feuillebois (19.89). All of these two-bubble papers 
noted the remarkable result that two equi-sized bubbles move with the same velocity 
as one bubble; Mevyappan et al. (1983) showed this numerically for the case of 
axisymmetric motion and Feuillebois proved it analytically, while Meyvappan & 
Subramanian and Anderson bowed that, for bubbles at any angle to the applied 
temperature gradient, the result was true to the order to which thev carried out their 
approximations. The study of electrophoretic phenomena, on the other hand, has 
important applications in areas concerned with the fractionation of mixtures of 
proteins or biological cell populations. Recall that electrophoresis has a venerable 
history, beginning with the work of Smoluchowski in the early part of this century: 
see Hunter (1981) for an extensive discussion of particle electrophoresis. But. for our 
purpose, the work of Reed & Morrison (1976) is especially noteworthy in that they 
showed that two identical spheres with thin double layers would move with the same 
velocity in an electric field as one sphere (the Smoluchowski velocity), irrespective of 
their separation and orientation in the electric field. 

In what follows, we shall extend the studies referred to above and explore their 
consequences for the averaging of thermocapillary and electrophoretic velocities in 
clouds of particles. The presentation is organized as follows. For the sake of 
simplicity, we first concentrate on the thermocapillary problem and then show how 
the electrophoresis problem can be addressed using the same mathematics. Thus, in 
§2 thermocapillary motion is discussed in detail and we show how particles in a cloud 
translate with the same velocity that one bubble would have if it were all alone. 
Thermocapillary motion of a space-filling cloud and the effect of backflow are then 
taken up in §3 where the renormalization ambiguity is set forth in detail. Next, the 
formal series for the many- body interactions is developed in §4. where we discover 
how to reduce the integral based on three- body interactions to one based on two- 
body interactions. Once the thermocapiilarv problem has been discussed, it is a 
simple matter to show how the analysis can be adapted for electrophoresis. This is 
done in §5 where we show that particles in a cloud of identical spheres move at the 
same velocity as an isolated sphere. In §6. we derive a surprising result which relates 
uniquely the average thermocapillary or electrophoretic velocity of a test particle to 
the effective conductivity of the suspension over the whole range of the particle 
concentration c. The paper concludes with a discussion in §7. 

2. Equations of motion 

Consider, first, a bubble in an infinite expanse of fluid of viscosity y and density 
p. in which a uniform, constant temperature gradient H far from the bubble is 



Particle migration in suspensions 99 

imposed. Let the surface tension of the interface between the bubble and the fluid be 
y. and assume that y decreases linearly with temperature T, and that the surface 
tension is high enough to keep the bubble spherical (of radius a). Then, as shown by 
\oung et al, (1959), the bubble moves with a velocity 


^vgb 0 


d yH 
d T'lp 


under the assumptions of low Reynolds number 


R — paL ygb//^ 


and small Marangoni number 

Ma 


al \-c,B = H &y q 2 
a dT'lpcc 


where a is the thermal diffusivity of the fluid. 

W e can use these results to non-dimensionalize the equations governing the 
temperature and velocity fields in this limit. Distances are scaled using the bubble 
radius a. the temperature is scaled using aH and the bubble centre is placed at the 
origin. We shall denote the dimensionless temperature gradient by h = H/H. The 
temperature field T is then determined by 


V 2 T = 0, 

together with the boundary conditions 


( 2 . 1 ) 


VT->h as |je| -*x, (2.2) 

and n-VT = 0 on |*| - 1 . ( 2 . 3 ) 

The velocity and pressure fields are scaled with the Young-Goldstein-Block velocity, 
giving 

— Vp + V 2 u = 0, ( 2 . 4 ) 

and \7 -h = 0, (2.5) 

together with the boundary conditions 

0 as (*| ->co. ( 2 . 6 ) 

un=Uh>n on |*| = 1 , ( 2 . 7 ) 

and o n-nn o n - 2V7 1 on|*| = l. (2.8) 


In (2.7) the unknown velocity of the bubble is written as Uh. but since we have scaled 
velocities with the known solution. V will equal 1 for the single-bubble problem. The 
boundary condition (2.8) has been simplified from its usual form using (2.3). The 
solution to these equations, for \x\ ^ 1 . is 


T = h-x-\h-V\xr, (2.9) 

and u = - \h- W, ( 2 . 10 ) 

where J is the Green function for the Stokes equations : 

J(x) = //|*| + x*/j*| 3 . 

The solution can be interpreted as follows. From (2.9) we see that the response of the 
bubble to the temperature gradient is equivalent to adding a thermal dipole at the 
bubble centre. Also, a bubble held fixed in a temperature gradient exerts a force on 
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the fluid (Subramanian 1985), but our problem requires the bubble to be force free, 
and therefore it translates with whatever velocity is needed to create a cancelling 
force. The difference between the two flow fields is a Stokes quadrupole, which is 
represented by (2.10). 

We now wish to consider the problem of n interacting bubbles, all of the same size, 
in an infinite expanse of fluid. We shall suppose that the centre of bubble p is at r p 
and let x v be a vector from the centre of that bubble to an arbitrary point x. so that 
x — r p + x p . Boundary conditions must now be applied on each surface |jc | = 1. We 
shall denote the velocity of bubble p by U p . The case of two bubbles (centres at r 1 
and r 2 ) was studied by Anderson (1985a) using the method of reflections, and he 
found that the velocities were given by 

Ui “ U 2 = h + 0( |x a -x 2 r 8 ). 

In other words, the velocity of one sphere was unaltered by the presence of a second 
to the indicated order of the approximation. We begin by extending Anderson’s 
analysis to three spheres and arrive at the equivalent result. In order to do so, we first 
rewrite some of his intermediate expressions in a more compact form. Anderson 
found that in addition to (2.9) and (2.10) he needed the velocity and temperature 
fields around a bubble in a quadratic temperature field, obtained from the ambient 
field T by evaluating WT at the point r p . In dimensional variables, we can write 

T= (VVT) p :(xx + $W\x\- 1 ), (2.11) 

u = ^r^ ( VV7’) p .-V( 1 +iV 2 )J. (2.12) 

As with the constant-temperature-gradient case, we can interpret these results using 
multipoles. A bubble responds to a linear temperature gradient by generating a 
quadrupole thermal disturbance and a stresslet (Stokes dipole) velocity disturbance. 
It should be noted that the stresslet flow field decays only as |*| -2 whereas the flow 
field (2.10) decays as |jc| - 3 , which is why it cannot be neglected in interaction 
calculations. 

We now consider three bubbles with centres at r v r 2 and r 3 . Starting with the 
bubble r x in the temperature gradient h . we find the disturbance temperature field at 
bubble r 2 from (2.9). Letting r 21 = r 2 — we find 

VT 21 = lh-(3r n rJ\r u \s-//\r n n 

This ambient temperature gradient induces a thermal dipole in sphere r 2 and gives 
it an additional velocity equal to VT 21 . However, the velocity field (2.10) around 
sphere r x creates an ambient velocity field at bubble r 2 equal to 

u 2 i = ~ih'(Zr 2l r 21 /\r 2l \ s ~//\r 2l \ 3 ). 

Thus the total change in the velocity of bubble r 2 is VT 2l + u 21 = 0. We now continue 
the method of reflections to sphere r 3 . 

The velocity of bubble r 3 will consist of contributions from the disturbances arising 
from the applied temperature gradient being reflected from bubble r x and from r 2 
together with the disturbances being reflected from bubble r x to r 2 and then on to r 3 . 
The last effect is the one we are interested in. We saw above that bubble r x induces 
a thermal dipole in bubble r 2 . This will produce a change in the ambient gradient near 
r 2 equal to 


iVT 21 -(3r 32 r3 2 /|r 32 ! 5 -//|r 32 l 3 ). 
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Bubble r 3 will therefore pick up a velocity equal to this. Also there is a relative 
velocity between bubble r t and the ambient flow at r 2 equal to VT 21 , and this will 
cause a velocity disturbance at r 3 . The ambient velocity at r 3 will be changed by 

-}vr 21 .(3r 32 r 32 /|r 3I P— //|r 32 |>), 

which cancels the temperature-gradient contribution. Finally we must consider the 
perturbations due to the Stokes dipole field given in (2.12). There are two 
contributions to this, of which the first comes from the second derivative of the 
temperature field (2.9) and the other from the rate-of-strain field obtained from the 
flow field (2.10). It has been shown by Anderson that these also cancel. 

We have seen then that two and three spheres move in an applied temperature 
gradient with velocities unchanged from the single-sphere case. We now wish to 
extend this to n spheres. Returning to (2.10), we can rewrite it using standard 
identities as 

u = \h- VV|*r i = V(|/r- V|*r j ). (2.13) 

From this, we see that if we write (2.9) for the temperature field as 

T = h'X + <fi, 

we can rewrite (2.10) as u = — V<}>. We now postulate that the last two expressions 


hold true for the flow around n spherical bubbles, and replace (2.1)-(2.8) with the 
equations 

T = h-x + <f>, (2.14) 

and « = -V0, (2.15) 

where = 0, (2.16) 

together with the boundary conditions 

V0-*O as |jc|h>-oo, (2.17) 

n'V(j) = —h-n on |* p | = 1. (2.18) 

The velocity field (2.15) will automatically satisfy the Stokes equations with the 
pressure equal to zero. In addition, the kinematic condition on the spheres’ surfaces, 
that sphere p moves with velocit}’ U p , is 

un — U p n = — hn on |x p | = 1, (2.19) 

meaning U v = h for all p. The boundar\ T condition left to verify is the stress 

boundary condition 

— n-(VV^)-f #tfi*(W0)*n = h + V<j> on |jc p | = 1 ' (2.20) 


To prove this we use an expansion of <p about r p , the centre of sphere p. To introduce 
this expansion we consider first the special case of just two spheres. The method of 
twin multipole expansions calculates (f> for this case by expressing it as a superposition 
of multipole expansions (Jeffrey 1973). The expression for (j> is valid everywhere 
outside the spheres, but it is expressed in more than one coordinate system (Jeffrey 
1973. equation 5.2). The coefficients in the expansion are found by transforming it 
entirely into the coordinate system of one sphere using addition theorems for 
spherical harmonics (Jeffrey 1973. equation 5.3) and applying boundary conditions. 
Here, although our aim is not to calculate coefficients, but rather to verify (2.20), we 
can still use the transformation idea extended to n spheres. Now <j> will be a 
superposition of n multipole series in decaying harmonics, with each series centred on 
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one of the n spheres. But each multipole series can be transferred to the centre of 
sphere p. where it will become a series of growing harmonics, just as can be seen in 
the two-sphere case. Thus <p can be expanded about r p in an expansion of growing 
and decaying harmonics valid in a spherical shell centred at r p whose inner surface 
is \x p \ = 1 and whose outer surface touches the surface of the nearest neighbour to r 
Let us denote a general decaying harmonic of order —n— 1 by . Then 
W 0-n-i "HI be a growing harmonic of order n. V\ e now write the expansion of 
$ about r p as 

f) — — ^•V|* p |~ 1 + £[(tt + l)|x p ! 2n + 1 + n ]^_ |j _ i< (2.21) 

* n 

In this form, (2.18) is automatically satisfied. Using the identity 

x p 'V(j>_ n _ x = ( — n— (2.22) 

together with the fact that x p = n on |jt p | = 1, one obtains, on |jc p | = 1, 

v <> = — \nh-n + \h + ^? i (2n+ l)(n + l)/i0_ n _ 1 + (2w + 1) 

n 

/i-(VV0) = lnh-n-lh + ^(2n + l)(n+l)(n~l)nf>_ n _ 1 ~(2n + l)Vf>_ n _ 1 . 

n 

and 

#i/»-(VV0)-n = 3nA-/i + 2(2w + l)(n + l)/i/i0_ n _ 1 , 

n 

whence one can easily verify that (2.20) holds. Thus the general result is proved. 
Incidentally, since <j> is harmonic, (2.15) satisfies the full Xavier-Stokes equations, 
with the pressure given by Bernouilli s equation, as long as the Marangoni number 
remains small. This exact solution to the thermocapillarv problem was also noted 
earlier for the single bubble case by Crespo & Manuel (1983) and Balasubramaniam 
& Chai (1987). 

3, The average velocity of a bubble in a cloud 

Since the velocity of a bubble is unaltered by the presence of other bubbles, it 
might seem that there is no averaging problem to perform. Thus the first 
renormalization strategy would be to suppose that U — U YCB . This is unsatisfactory, 
however, as can be shown by a physical argument-. If the cloud of bubbles fills a 
closed container, conservation of volume (the fluid is incompressible) show's that the 
velocity must change, for if w r e take a surface through the suspensions perpendicular 
to the applied temperature gradient, the rate at which volume is displaced by the 
bubbles equals cU YGB , w here c is the volume concentration of the bubbles. Since the 
net flux across the surface must be zero, there must be a net backflow of fluid 
reducing the velocity of the bubbles by a factor (1-c). 

A second argument in support of renormalization comes from noting that the exact 
cancellation of the two velocity perturbations only occurs for non-conducting 
bubbles of the same size, because, as Anderson (1985a) showed, if the bubbles (or 
drops) have any conductivity, or if they differ in size, there is a familiar 1 /|jc| 3 non- 
convergent interaction between them. Thus, for the more general conducting-drop 
problem, we must renormalize, and then if w r e take the limit of zero conductivity, we 
are led to a renormalized integral, not the result U = U yGB . However, although these 
arguments help to settle the question of whether or not renormalization is needed. w r e 
must still decide on the renormalization procedure. Recall that from the overall 
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specification of the problem we have constraints on u{r). the velocity at the origin due 
to a bubble at r, and VT(r), the temperature gradient at the origin. They are 

juP(r)dr = 0 and j(VT-H)P(r)dr = 0. 

Both of these are non-convergent integrals that cannot be evaluated From the 
results above, however, it can be shown that the two integrands have the same 
dependence on r when the bubbles do not overlap the origin (r > 1 ) but thev are 
different when the bubbles overlap the origin. Therefore the integral 

J ([ U(r) - t/ YGB ] P(r)0) - [A, u 4- A 2 ( VT- H)] P(r )) dr 

will be convergent for any value of A x , provided A 1 = A 2 , since both terms in brackets 
will then be zero for r > 2. The expression evaluates to -f \ x cH and therefore the 
mean \ elocity of the particles in the suspension is 

( 1 f Aj c) C^ygb 

if the particles are non-conducting. Anderson’s result is obtained by setting A 1 = 1. 
In fact, even when the particles are conducting or have different ‘sizes, there still 
exists an indeterminacy, because A x and A 2 must satisfy onlv one equation, namely 
that the non-convergent terms in the two sets of brackets cancel. Analogous 
considerations hold for particles undergoing electrophoresis (Anderson 1986). 

Chen & Acrivos (1978) discovered a similar indeterminate situation, by finding 
several possible convergent two-body integrals. As here, the integral could 
apparently be put in a convergent form without renormalization and two 
renormalized integrals were found. Since these authors were able to resolve the 
indeterminacy by considering three-body 7 interactions, w r here onlv one choice gave a 
convergent integral, we now extend the group expansion approach of Jeffrey (1974) 
to include the extra field present, in order to investigate the properties of the 
expansion and resolve the indeterminacv. 


4. Formal series for three-body interactions 

In order to write down the integral for three-body interactions, we must add to the 
notation in Jeffrey (1974) the fact that the velocity of a bubble in the laboratory 
frame now depends upon two applied fields, viz. the ambient temperature gradient 
and the ambient velocity. We begin by considering a test bubble at the origin 
surrounded by n others, the whole cloud being immersed in a temperature gradient 
Ff and an ambient velocity 7 u. At the end of the calculation, we shall set u to zero, but 
we have seen in the method-of-reflections solution above that we must -consider 
changes in the ambient flow. We shall write n \H\U) as the general notation for 
the velocity of the test bubble, <€ n being the set of position vectors r p of the other 
bubbles. We define incremental velocities in the manner of Jeffrey (1974). Thus 
U 0 {-H:u) is the velocity of a single sphere in an ambient temperature gradient H 
and a velocity 7 field u. With suitable non-dimensionalization we know 

U 0 (;H-O) = H+u. 

We next consider the test bubble at the origin with a second bubble (a first 
neighbour) at r r The velocity of the test bubble is now Uir^.H-.u). From this we 
define the increment U 1 due to the second sphere as 

t/ 1 (r 1 ;tf;w) = U{r , : H; u) - U 0 ( : H: a). 
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Xow since two bubbles respond to an ambient velocity a bv each increasing its 
velocity by the same speed u. it follows that £/, is independent of a. Therefore 

= U 1 (r 1 ; R;0) = V(r, ;R;0)~ U a (-R;0). 

In addition we have seen that two bubbles move at the same speed as one, so in 
this problem U x = 0. but in general this would not be so. 

The continuation to a third bubble (neighbours at r x and r 2 ) uses the obvious 
notation 

t/ 2 (r, r.. Rm) = t/(r, r 2 ; H, a ) - lf,( r t ,R,a)~ \R\a)~ U 0 (;R; a). 

Again we see that U, is independent of u and, for our bubble problem is actually 
identically zero. 

The first group expansion is the generalization of the non -renormalized integral It 
consists simply of averaging the velocity increments just defined : 

U= U 0 + ju i (r l -,R-a)P(r l \0)dr l + jju i P{r l r 2 \ 0 )dr 1 dr i + .... 

Since U x = U 2 = 0, etc. we again arrive at the unphysical result U = U = U Y 
This is different from the previously known case of Chen & Acrivos (1978) in which 
the expansion broke down visibly at the three-bodv term bv giving a non- convergent 
integral. " 6 

The second group expansion corresponds to generalizing the renormalized integral. 
To introduce the notation required to write this down, we shall ignore for the 
moment the fact that U x = U 2 = 0. It is not true in the more general case of drops 
with non-zero conductivity, and we need to write down the general asvmptotic 
relations which hold when the spheres are far apart in order to obtain the second 
expansion. The latter is obtained by rewriting the method-of-reflections results 
above in a general notation. To do this we need a notation for fields near the origin 
produced by distant spheres. Consider first a sphere by itself at r v The temperature 
gradient at the origin is H(r x :H), meaning the field produced by a sphere at r, 
immersed in an ambient gradient H. \\hen there are no spheres present, the field 
at the origin is just the applied field H. so we can define an increment in H to be 

= H ( r J H )r H This is the field iV(l/|*|) given in (2.9). The velocity field 
at the origin is \u{r x ;H:u). meaning a sphere at r x moving because of the field i?and 
ambient flow u. Again there will be a velocity increment: 

u x {r x :H-u) = u(r x :H,u)~u. 

This will be independent of u, so 

u x {r x \ H . u) = u x {r x ;H;0). 

The two- body integral of the second group expansion is now the obvious 
generalization of Jeffrey (1974): 

J{ U x (r x ; H ; u)P(r x \ 0) - U 0 ( ] H x (r x ; H); u x (r x ; H, 0)) P(r x )} d r x . 

This is the integral evaluated by Anderson. Since U x is zero, it reduces to a one- 
particle integral and gives the result 

C=(l-ic)C yGB . 

We now need to consider two spheres far from the test sphere. We can extend our 
notation in the obvious way ._H 2 {r x r 2 :H) is the increment at the origin due to spheres 
at r x and r 2 , while u 2 (r x r 2 \H\u) is the increment in the fluid velocity at the origin. 



105 


Particle migration in suspensions 

W e shall not repeat the elaborate asymptotic arguments needed to establish the 
three-body term in the second (renormalized) expansion, but write it down as the 
obvious extension of Jeffrey (1974). It is 

I ( U 2 (r 1 r 2 : H: u) P(r x r 2 \ 0) - £/, (r, u x {r 2 :H, 0)) P(r x \ 0) P(r 2 ) 

- U 0 ( : H 2 {r 1 r 2 \H) ; u 2 (r 1 r 2 ; H ; 0)) P{r 1 r 2 ) 

+ U 0 ( ; H 1 (r 1 ; H 1 (r 2 :H))\ u^{r l ;H 1 {r 2 -,H)\ 0)) Pf/q) P(r 2 )} dr l d r 2 . 

We can now integrate over r i and r 2 and obtain the 0{c 2 ) term. First we simplify. We 
know that U 2 — U 1 = 0. So we need consider only the last two terms. 

Since u = — V0 outside the particles, the third term U 0 ( ; H 2 ; u 2 ) will be zero unless 
one of the particles overlaps the origin, and in that case u 2 - 0 because the velocity 
of two bubbles equals the velocity of one particle. Therefore the third term is non- 
zero only if IrJ ^ 1, when we have 


- U Q (\H 2 \u 2 )P(r^ r 2 )= -H 2 {r l r 2 :H)P(r 2 1 r^P^). 

Now we perform the integration by fixing r 2 — r l = s and integrating first over r 1 
followed by integration over s. The integral over r l is just the volume integral of H 2 
over the spherical bubble, which equals the thermal dipole strength of the bubble S 1 . 
The integral becomes 


H 2 {r l r 2 :H)P(r 1 r 2 )dr 1 dr 2 - S 1 {r 2 :H)P{r 2 \r 1 )dr, 


Finally the last term in the three-body integral can be treated the same way to 
become r 


S 0 (:H 1 (r 2 :H))P{r 2 )dr 2 


Hence the three-body velocity problem reduces to the two-bodv heat conduction 
problem. In §6. we shall show that this is part of a more general result. 


5. The electrophoresis problem 

To show that much of the analysis and all of the conclusions derived thus far 
for thermocapillarity apply equally well to charged spheres moving through an 
electrolyte in an electric field we first recall that a single, non-conducting sphere 
moves in a uniform electric field E at the Smoluchowski velocity 



as long as the double layer is thin (cf. Hunter 1981). Here e is the dielectric constant 
of the suspending fluid. e 0 is the permittivity of free space, £ is the electrostatic 
potential at the surface of the particle (the zeta potential), and p is as before the 
viscosity of the suspending electrolyte. Equation (o.l) applies when the particle 
motion is slow and the diffuse charge thickness (the Debye length a: -1 ) is small 
compared with the particle radius a. Viewed on the lengthscale of the particle, fluid 
appears to slip past the surface at a speed known as the electro-osmotic slip velocity. 
However, when the details within the thin region adjacent to the surface are resolved 
by. for example, singular perturbation methods, the no-slip condition holds and the 
velocity in the diffuse region joins smoothly to that outside. Inside the diffuse layer 
the Stokes equations must be modified by adding a term to account for the 
electrostatic body force on the fluid due to the action of the field on the space charge. 
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For thin double layers and non-conducting particles, however, we can avoid 
analysing events within the diffuse layer in a consistent fashion by replacing the no- 
slip condition bv „ 

rr ee 0 £_ . 

U = u—e ; -Vt (5.2) 

and using the Stokes equations without the extra body force. Here U denotes the 
particle velocity, u is the fluid velocity, and - i s the local value of the electric field 
(Anderson 19856; Russel. Saville & Schowalter 1989). 

The electric field is governed by solutions of Laplaces equation with the normal 
gradient equal to zero at a particle surface (non-conducting particles). Inside each 
particle. Laplace’s equation applies and the two fields are continuous at the interface. 
The potential gradient is constant, i.e. -Vij/^Eh. far from the particle, while the 
velocity and pressure fields follow from solutions of Stokes equations with the 
velocity vanishing far from th^ particle. The boundary condition at the particle 
surface, (5.2), shows that the normal velocity of fluid and particle match at the 
surface. 

Now, the electrophoresis problem can be restated in a form analogous to that used 
in §2 for thermocapillary motion. First, scale the potential on ( — aE) and recover 
(2.1 )— (2.3) to describe the dimensionless potential outside the sphere. Next scale 
the velocity on the Smoluchowski velocity. ee 0 £E/ju, to obtain (2.4) and (2.5). The 
boundary condition (2.6) still holds but in place of (2.7) and (2.8) we use the 
dimensionless form of (5.2), namely 

U=u + VT. ( 5 . 3 ) 

It follows immediately that the dimensionless potential and velocitv are given bv 
(2.9) and (2.10). 

With these results for a single sphere in hand we can investigate the electrophoresis 
of two identical spheres. Here, as long as the separation exceeds the Debye thickness, 
the spheres move at the Smoluchowski velocity irrespective of their separation and 
orientation (Reed & Morrison 1976). This result also holds for n identical spheres 
(Zukoski 1984). as can be seen using the arguments given at the end of §2 if (2.8) 
is replaced bv 

U p h = u + VT on |* p | = 1. (5.4) 

It follows that L p = 1 and the expansion for $ in the harmonic series noted earlier 
holds. Accordingly, we have the somewhat surprising result that each particle in a 
cloud of spheres undergoing electrophoresis moves at the same velocity. 


6. Equivalence of the migration velocity and the effective conductivity 

It was shown in §4 that tiro integral giving the three-body correction to the 
thermocapillary or the electrophoretic migration velocity can be reduced to the 
integral giving the two-body correction to the effective conductivity of the 
suspension. Here we shall derive a remarkable result that the migration velocity can 
be obtained simply from the effective thermal conductivity of the suspension of non- 
conducting particles for all particle concentrations. Specifically, taking the gradient 
of (2.14). we observe that VT = h + T/( j ) (6 j) 

and averaging this over the volume of the suspension V, we obtain 



VTdF 


= h. 


( 6 . 2 ) 
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Splitting the domain of integration into F f . the part of V occupied by fluid, and V . 
the part occupied by particles, we arrive at 

h= (l-c)h+yj VfSdr + y;^ VTdF. (6.3) 

In the notation of Jeffrey (1973). the last integral is simply equal to — cS . the average 
thermal dipole strength of the particles. 

Before we can a\ erage (2.15) we must remember that this is for particles in an 
infinite fluid, whereas in a container a backflow V B will exist which must be included 
in the equation. Thus (2.15) is replaced by 

u{x) — — V$+U B , when x is outside a particle, (6.4a) 

u{x) — V— h+ U B , when x is inside a particle. (6.46) 

Averaging this equation, we obtain 

e = 0 = --ij V0dT’ + (1 — r) U B + cU. (6.5) 

l sing (6.3) to eliminate the unknown integral from (6.5) and then using (6.46) to 
eliminate U B from (6.5). we obtain 

17= h + cS = K*-h, 

where K* is the effective conductivity of the medium. It should be noted that this 
result holds for any particle concentration, and allows us to use not only the existing 
direct calculations of K* but also bounds on the effective conductivity, such as the 
Hashin-Shtrikman bounds (Hashin & Shtrikman 1962) or more recently developed 
expressions (Torquato 1987). For example, using the results for an isotropic well- 
mixed suspension given in Jeffrey (1973). we obtain 

U = (1 — §r + 0.59c 2 + ...) h 

for the same suspension. 


7. Discussion 

The problems of electrophoresis and thermocapmarity provide an interesting 
challenge to the established methods for calculating the average properties of 
suspensions. Although it is possible to take refuge in the more general cases of 
conducting particles or particles of different sizes to distinguish between the 
alternatives here (in conjunction with the three-body problem), this is an unsatisfying 
solution from a general point of view since there is no guarantee that such 
generalizations will always be possible. Moreover, special cases are usually chosen for 
study because they are simpler than the general problem, so it is disappointing to 
have to return to a general situation in order to address the special case. One 
argument used successfully in the past to choose between alternative schemes 
consisted of appealing to higher approximations. Thus, conflicting convergent 
integral expressions for the two-body interaction were resolved bv extending the 
analysis to three-body schemes. It seemed from the case of Chen & Acrivos (1978) 
that the second group expansion of Jeffrey (1974) would always be able to distinguish 
between different possible treatments of the two-body problem by jumping to the 
three- body problem (another example of generalizing the problem). The particle 
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problems studied here, however, show that this is not so. because the entire series for 
the first and second group expansions can be searched without finding an integral 
that converges in one case but not the other. Thus if we do not wish to leave the 
confines of the non-conducting particle problem, there is only the argument based on 
conservation of mass to choose between the mathematically possible alternatives. 

Experimental tests of the theory for the thermocapiilary or electrophoretic 
velocity of particles in a suspension are extremely difficult. To begin with, problems 
arise because of buoyancy and polydispersity. In addition most systems become 
opaque well before the particle number density is large enough to influence the 
average translational velocity and this makes it hard to measure the speed of 
individual particles. Nevertheless, Zukoski & Saville (1987, 1989) were able to 
prepare neutrally buoyant, transparent suspensions by removing the hemoglobin 
from red blood cells. A variety of particle shapes could be formed. Then tracer 
particles with the same electrokinetic properties as the red-blood-cell ghosts could be 
tracked microscopically over a wide range of particle concentrations. Zukoski & 
Saville measured the electrophoretic mobility of individual particles and the effective 
conductivity of the suspension over volume fractions ranging between 0 and 0.8. The 
effective conductivity was found to follow the Maxwell relation for slightly 
conducting particles while the mobility varied as (1— c) up to the highest volume 
fractions. The reasons for the difference between their results and the theory 
developed here are unclear. It may be due to the slight conductivity of the particles, 
but further work is obviously necessary. 
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Postscript 

The references in the text acknowledge formally the influence Professor Batchelor 
had on this work, but, less formally, the interplay between physical and mathematical 
arguments discussed in this paper brings to the- mind of one of us (D.J.J.) an 
exchange that took place during a seminar at Cambridge in the mid 1970s. 
speaker (justifying elaborate mathematical argument) : The trouble with a physical 
argument is that you may not get all the terms. 

GKB : No. that is the trouble with a bad physical argument. 

several voices: How do you tell a good physical argument from a bad one? 

GBK : That's easy: you think. 
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